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Blowout  is  due  to  severe  overheating,  where  the  inside  of  the  pad  appears  to  vaporize. 

The  properties  of  the  pad  which  cause  temperature  related  problems  are  low  thermal 
conductivity,  high  compliance,  and  high  internal  damping.  The  second  and  third  propert 
are  important  to  the  basic  functions  of  the  pad,  so  care  must  be  taken  that  they  are 
maintained  to  some  extent  in  any  redesign.  The  approach  taken  here  is  to  modify  the  ba: 
rubber  pad  design  by  adding  reinforcement.  Metal  wire  reinforcement  %was  chosen  because 
it  will  help  directly  by  strengthening,  and  indirectly  by  increasing  the  conduction  of 
heat  to  the  pad  surface.  To  maintain  compliance  and  to  add  strength  in  the  most  critic, 
direction,  the  mode  of  reinforcement  chosen  was  uniaxial,  in  the  direction  of  travel  of 
the  road  wheels.  The  recommended  baseline  design  was  a  standard  T142  configuration  witl 
brass  or  bronze  coated  steel  wire  of  32  mil  diameter  and  a  volume  fraction  of  0.013, 

Analytical  work  included  a' thermal  model  of  a  track  pad,  which  quantified  the  reduction 
in  operating  temperature  due  to  reinforcement.  It  was  also  used  to  model  the  curing 
process.  A  fracture  mechanics  based  fatigue  model  was  constructed,  to  model  debonding 
of  the  reinforcement  as  a  function  of  cyclic  loading.  This  showed  that  the  baseline  des 
should  maintain  a  bond  for  a  reasonable  mileage,  and  showed  what  design  changes  would  b< 
needed  to  improve  performance.  Adhesion  tests  were  conducted  for  a  range  of  rubber 
compounds  and  wire  types,  and  likely  combinations  with  high  adhesion  were  identified. 
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1.0  INTRODUCTION 


Current  tank-track  designs  use  rubber  pads  (track  pads)  to  enhance  traction, 
vibration  Isolation,  and  noise  reduction,  and  to  limit  damage  to  paved 
surfaces.  The  pads  can  be  made  Integral  with  the  track  (as  In  the  M1/T156 
system),  or  bonded  to  steel  plates  which  are  then  bolted  to  the  track  (as  In 
the  M60/T142  system).  Typical  track  pad  mileage  Is  far  below  the  mileage 
between  major  overhauls.  Since  pad  replacement  Is  expensive  and  time 
consuming,  ways  of  Improving  the  design  to  Increase  mileage  while  maintaining 
the  functions  stated  above  are  needed. 

Failure  mechanisms  for  track  pads  Include  abrasion,  cutting,  chunking,  and 
blowout.  All  of  these  are  related  to  hysteresis  heating  of  the  rubber,  which 
Is  due  to  the  cyclic  loading  of  the  road  wheels  on  the  track.  Rubber  loses 
tensile  and  tear  strength  rapidly  as  temperature  Is  Increased.  Cutting  and 
chunking  occur  oh  rough  terrain,  when  pads  are  loaded  unevenly  by  the  sharp 
edges  of  rocks.  Blowout  Is  due  to  severe  overheating,  In  which  the  Inside  of 
the  pad  appears  to  vaporize. 

Pad  properties  that  cause  temperature  related  problems  are  low  thermal 
conductivity,  high  compliance,  and  high  Internal  damping.  The  second  and 
third  properties  also  are  Important  to  the  basic  functions  of  the  pad,  so  care 
must  be  taken  that  they  are  maintained  to  some  extent  In  any  redesign.  The 
approach  taken  here  Is  to  modify  the  basic  rubber  pad  design  by  adding 
reinforcement.  Metal  wire  reinforcement  was  chosen  because  It  will  help 
directly  by  strengthening,  and  Indirectly  by  Increasing  the  conduction  of  heat 
to  the  pad  surface.  To  maintain  compliance  and  to  add  strength  In  the" most 
critical  direction,  the  mode  of  reinforcement  chosen  was  uniaxial,  In  the 
direction  of  travel  of  the  road  wheels.  The  recommended  baseline  design  was  a 
standard  T142  configuration  with  brass  or  bronze  coated  steel  wire  of  32  mil 
(0.8  mm)  diameter  and  a  volume  fraction  of  0.013. 

Analytical  work  Included  a  thermal  model  of  a  track  pad,  which  quantified  the 
reduction  In  operating  temperature  due  to  reinforcement.  It  was  also  used  to 
model  the  curing  process.  A  fracture  mechanics-based  fatigue  model  was 
constructed  to  model  debonding  of  the  reinforcement  as  a  function  of  cyclic 
loading.  This  showed  that  the  baseline  design  should  maintain  a  bond  for  a 
reasonable  mileage,  and  showed  what  design  changes  would  be  needed  to  Improve 
performance.  Adhesion  tests  were  conducted  for  a  range  of  rubber  compounds 
and  wire  types,  and  likely  combinations  with  high  adhesion  were  Identified. 

2.0  OBJECTIVE 

4 

The  general  aim  of  this  program  was  to  extend  the  service  life  of  tank  track 
pads.  The  specific  goals  were  to  1)  to  develop  a. baseline  design  for  a  fiber 
reinforced  T142  type  track  pad,  2)  to  present  reasonable  mathematical  and 
physical  bases  for  the  design,  considering  thermal  and  adhesion  behavior,  3) 
to  carry  out  a  series  of  adhesion  tests,  to  Identify  promising  combinations  of 
rubber  compounds  and  reinforcement,  and  4)  to  model  the  molding  process,  so 
that  ways  to  achieve  more  uniform  properties  by  reducing  over  and  under  curing 
could  be  presented. 
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3.0  CONCLUSIONS 


An  analytical  basis  has  been  presented  to  support  Improvements  In  track  pad 
life  via  steel  wire  reinforcement.  It  Included  a  simple  1-0  thermal  analysis 
which  gave  results  consistent  with  field  data  from  the  literature,  a  fracture 
mechanics  based  study  of  the  Influence  of  temperature  and  reinforcement 
parameters  on  tearing  and  wire  debonding,  a  wire  pullout  test  program,  and  2-D 
numerical  models  for  hysteresis  heating  and  curing.  Stress  In  the  rubber  can 
be  reduced  by  orienting  the  wire  In  the  direction  of  highest  tensile  stress, 
which  Is  the  direction  of  road  wheel  travel.  One-dimensional  reinforcement 
allows  the  pad  to  retain  most  of  Its  original  compliance.  Temperature  In  the 
pad  Is  reduced  through  stiffening,  and  through  Increased  thermal  conductivity 
In  the  wire  direction.  Additional  temperature  reduction  can  be  achieved  by 
Improved  surface  heat  transfer.  By  Itself  (l.e..  Ignoring  reinforcement 
effects),  the  resulting  temperature  drop  Is  large  enough  to  reduce  the  rate  of 
tear  growth  In  the  rubber  by  a  factor  on  the  order  of  1/2. 

The  Issues  of  wire  strength  and  bond  strength  have  been  addressed.  Wire 
stress  at  a  reasonable  volume  fraction  requires  use  of  steel  wire.  (Brass, 
aluminum,  and  copper  are  not  strong  enough.)  The  maintenance  of  a  bond 
between  the  wire  and  the  rubber  was  recognized  as  a  critical  Issue,  with 
failure  postulated  as  a  large  debond  length.  Bond  Integrity  under  cyclic 
loading  has  been  studied  using  a  fracture  mechanics  based  fatigue  model,  which 
was  correlated  to  the  results  of  single  cycle  wire  pullout  tests.  A  method 
has  been  presented  for  estimating  the  Improvement  potential  of  changes  In  wire 
size  and  volume  fraction.  A  test  method  from  the  literature  was  adapted  to 
the  nonlinear  behavior  exhibited  by  the  rubber.  Wire  pullout  tests  were 
conducted  as  part  of  this  program,  using  several  rubber  types,  and  several 
types  of  adhesive  and  metallic  coatings  on  the  wires.  The  strongest  bond  by 
far  was  achieved  with  a  high  sulfur,  natural  rubber  compound  and  brass  coated 
steel  wire.  Good  bonds  were  achieved  with  an  EV  cure  rubber,  using  adhesive 
on  a  copper  coated  steel  wire  and  on  a  brass  coated  steel  wire.  The  process 
of  applying  the  adhesive  required  on  etch  for  the  copper  surface,  but  not  for 
the  brass  surface.  In  any  case,  use  of  an  adhesive  proved  to  be  very 
cumbersome.  A  good  bond  was  also  achieved  using  a  TACOM  standard  compound 
(14A)  with  bronze  coated  steel  wire,  and  no  adhesive.  Because  adhesive  was 
not  needed,  and  because  the  rubber  compound  had  already  been  shown  to  perform 
well  In  track  pads,  this  combination  Is  recommended  for  testing. 

Numerical  studies  were  carried  out  with  a  two-dimensional  finite  difference 
thermal  analysis  code  written  specifically  for  this  project.  The  code  models 
hysteresis  heating  and  the  curing  process  for  track  pad  sized  rectangular 
blocks.  The  results  for  hysteresis  heating  showed  that  the  Increase  In 
thermal  conductivity  due  to  steel  wire  reinforcement  can  lead  to  moderate 
decreases  In  pad  temperature.  Stiffening  effects  can  be  Important  because 
they  directly  reduce  the  heat  Input  from  cyclic  loading.  Improvements  In 
surface  heat  transfer  can  lead  to  temperature  reductions  on  the  same  order  as 
those  due  to  Increased  thermal  conductivity.  For  a  baseline  case  with  steel 
wire  reinforcement,  a  numerical  model  of  cyclic  loading  showed  the  following 
reductions  In  temperature  rise  from  the  following  phenomena:’ 
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stiffening  25% 

enhanced  thermal  conductivity  15% 

enhanced  surface  heat  transfer  15% 

The  combined  effect  Is  a  46%  reduction  In  temperature  rise,  which  could  have  a 
significant  effect  In  reducing  tearing  and  blowout  failures.  Reduced 
temperature  also  helps  maintain  bond  strength. 

Curing  models  showed  that  the  Improved  thermal  conductivity  due  to  wire 
reinforcement  can  Improve  the  uniformity  of  the  cure  In  a  thick  rubber  track 
pad.  Other  means  of  Improving  cure  uniformity  were  Investigated.  Preheating 
and  post-heating  at  temperatures  below  the  mold  temperature  (by  30  to  40°C) 
were  shown  to  hold  promise,  and  work  by  Inducing  higher  temperature  for  longer 
time  at  the  pad  center. 


4.0  RECOMMENDATIONS 

4.1  Wire  Reinforced  Track  Pads 

Fabricate  wire  reinforced  track  pads  as  per  Figure  4-1.  Bronze  coated  wire 
(hose  wire)  and  TACOM  baseline  triblend  (compound  14A  MIL-T-1 1 491 C  AT)  are 
recommended.  Consider  slotted  sides  for  Improved  surface  heat  transfer. 
Baseline  wire  size  Is  32  mil  (0.8mm)  diameter,  and  baseline  volume  fraction  Is 
0.013.  Wire  spacing  Is  0.25  Inch  (6.4mm). 

4.2  Cure  Cycles 

If  nonuniformity  of  cure  Is  a  problem,  consider  revising  the  cure  cycle  by 
subjecting  the  rubber  to  an  appropriate  preheat  and/or  post-heat,  or  a  longer 
mold  time  at  a  lower  temperature.  Explore  the  cost  saving  possibilities  of 
reducing  the  required  mold  time  by  adding  a  post-heat,  out  of  the  mold,  at  a 
lower  temperature. 

5.0  DISCUSSION 

5.1  Background  Data 

Typical  mileage  for  M60/T142  track  pads,  as  reported  by  TACOM,  Is: 

2,000  miles  on  pavement 

900  miles  on  gravel 

250  miles  on  rocky  cross-country  terrain. 

In  contrast,  the  distance  run  between  major  overhauls  In  6,000  miles.  Costs 

as  reported  several  years  ago  were  $5  per  mile  for  fuel,  and  $27  per  mile  for 
track.  Clearly,  If  the  time  and  material  Involved  In  replacing  worn  track 
pads  could  be  reduced,  significant  savings  could  be  achieved.  A  goal  for  the 
long  term  would  be  to  Increase  pad  mileage  to  the  mileage  between  major 
overhauls . 


Measurements  of  temperature  rise  as  a  function  of  time  In  M60/T142  track  pads 
for  various  surfaces  have  been  reported1.  Approximate  values, 
extrapolated  to  steady  state  (more  than  two  hours  running  time),  are  given 
below: 


Surface 

Avg  Speed  (MPH) 

AT( Interior) ,°C 

AT(Surface) ,°C 

Paved 

20 

70 

50 

Gravel 

20 

90 

- 

Cross  Country 

16.7 

105 

- 

Under  desert  conditions,  the  Initial  temperature  of  the  pads  can  be  as  high  as 
50°C ,  due  to  the  high  ground  temperature.  Thus,  the  pad  Internal  temperature 
can  be  as  high  as  155°C  ( 3 1 0° F ) .  For  changes  In  vehicle  velocity,  the  rate  of 
hysteresis  heating,  and  thus  the  temperature  rise,  will  be  In  direct 
proportion  to  the  velocity. 

The  effect  of  temperature  on  tensile  strength  of  rubber  Is  shown  In  Figure 
5-12.  Tear  strength  and  tensile  strength  are  generally  closely  related. 

In  Figure  5-1,  the  material  that  most  closely  resembles  track  pad  material  Is 
the  Styrene-Butadiene  Rubber  (SBR)  compound.  Strength  Is  still  falling 
rapidly  as  temperature  rises  at  140°C. 

The  rate  of  tear  growth  as  a  function  of  temperature  under  given  loading 
conditions  may  be  a  good  way  to  judge  the  effect  of  temperature  on  track  pad 
mileage.  Standard  extrapolation  techniques  show  a  highly  nonlinear  positive 
relationship  between  rate  of  tear  growth  and  temperature.  Moderate  reductions 
In  temperature  could  therefore  result  In  large  Increases  In  pad  mileage. 

In  a  test  described  to  us  by  TACOM,  a  tank  was  run  over  a  glass  block  on  which 
a  metal  bar  was  placed  (to  model  a  rock).  As  the  road  wheels  rolled  over  the 
pad,  the  event  was  photographed  through  the  glass  plate.  This  record  clearly 
showed  high  tensile  strains  on  the  surface  of  the  pad,  In  the  direction  of 
travel  of  the  road  wheel.  The  deformations  Imply  a  large  tensile  stress  In 
the  Indicated  strain  direction.  This  helped  In  the  choice  of  a  direction  for 
the  reinforcement.  See  Figure  5-2. 

5.2  Analysis 

5.2.1  Influence  of  Temperature  on  Rate  of  Tearing 

Figure  5-1  shows  tensile  strength  as  a  function  of  temperature.  Tensile 
strength  Is  equivalent  to  the  fatigue  strength  with  one  cycle  to  failure.  The 
failure  mode  of  track  pads  Is  closer  to  a  many  cycle  fatigue  mode,  with  tears 
growing  a  small  Increment  with  each  cycle  of  loading.  A  way  of  evaluating  the 
effect  of  temperature  on  this  process  Is  a  WIlllams-landel-Ferry  (WLF) 
relation3,  fit  to  data  for  rates  of  tear  growth  at  temperatures  of  20°C 
and  70°C.  The  relation  can  be  expressed  as: 

R  =  rate  of  tear 
T  =  temperature 

Tg  =  glass  transition  temperature  (-63°C  for  SBR) 
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20  60  100  140 

(  .i  Temperoturt,  *C 

Fig.  3.3  Tensile  strength  v.  temperature 
•, Curve  1  Natural  rubber  gum  compound 

Curve  2  Natural  rubber  containing  50  parts  lampblack  per  100  of  rubber 
Curve  3  SBR  containing  50  parts  HAF  black  per  100  of  rubber 


Figure  J5-1.  Strength  Vs  Temperature  (Ref.  2) 


(1) 


RID. 

R(T0) 


a<y 

a(T) 


loglo[a(T)]  =  -8.86(T-Ts) 

101 .6+T-Ts 


T  =  T  +  20°C 
.  s  9 

Using  this  formulation,  the  effect  of  a  change  In  temperature  on  the  rate  of 
tear  growth  for  a  given  loading  can  be  estimated.  Figure  5-3  was  generated 
using  T0  =  160°C.  Reducing  the  temperature  from  160°C  to  120°C  results  In 
the  rate  of  tear  growth  dropping  by  64%.  These  temperatures  are  outside  the 
range  of  the  data.  Nevertheless,  the  trend  offers  some  useful  Insights. 

5.2.2  1-D  Thermal  Analysis  of  Track  Pad 

Much  can  be  learned  from  reducing  the  problem  to  Its  simplest  analytical  form, 
and  examining  a  closed  form  solution.  The  simplest  model  of  the  track  pad  Is 
one-dimensional,  through  the  thickness,  with  uniform  heat  generation  due  to 
hysteresis,  convective  boundary  conditions,  and  steady  state4.  A  summary 
Is  as  follows: 

W  =  heat  Input  rate  per  unit  volume  from  hysteresis 
T  =  temperature,  T$  *  T( surface),  Tj.  =  T(free  stream) 
k  =  thermal  conductivity 
h  =  surface  heat  transfer  coefficient 
Y  =  direction  of  heat  flux,  boundaries  at  Y  =  0,  L 
p  =  density 
c  =  specific  heat 
t  =  time 


pcaT  = 
at 


0  =  kd*T  +  W 


2 

dY 


-kdT  =  h(T  -T-)  at  Y  =  0,L 
dY  S  r 


U) 
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Choose  Fiber  Orientation  Based  on  Tensile  Stress,  Heat  Flow 


FIBERS  ARE;  CONTINUOUS  AND  UNIDIRECTIONAL 

Figure  5-2.  Orientation  of  fiber  reinforcement. 
Fibers  are  continuous  and  unidirectional 
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RELATIVE  TEAR  GROWTH  RATE  R/R(0) 


Figure  5-3.  Relative  tear  growth  rate  in  SBR  as  a  function  of 
temperature  (WLF  formulation) 
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Integration  yields 


T  =  WfY(L-Y)  ♦  L_,  +  Tf 

1  2k  2hJ 

AT  =  T  -  Tff  ATmax  =  AT( L/2)  =  ( WL/2) [ ( L/4k)  ♦  (1/h)] 

AT( Y=0,L)  =  WL 

2h  (3) 

Clearly,  Internal  temperature  will  be  reduced  If  W  Is  reduced,  or  If  k  or  h 
are  Increased. 

Temperatures  calculated  from  the  baseline  formulation  are  compared  below  to 
actual  measurements1*5.  Field  track  pad  temperature  measurements  were 
extrapolated  approximately  to  steady  state  for  comparison  purposes.  (The 
references  also  reported  numerical  predictions  made  using  a  detailed  finite 
element  model.)  The  relevant  Inputs  for  the  1-0  model  are  as  follows: 

5  _3  _1 

W  =  3. *10  erg  cm  s 
k  =  3.17*10*  erg  cnfV^C"1 

4  _2  _1  _ 1 

h  =  1 .7*10  erg  cm  s  °C 
L  =  2.125  Inches  =  5.4  cm, 
measurements  at  Y  =  0,  1.6  cm 

Using  the  1-0  equations, 

AT( Y=0)  =  48°C 

AT( Y=1 .6  cm)  =  76°C  (4) 

AT( Y=2 . 7  cm)  =  82°C 

Extrapolating  from  time-temperature  data  presented  In  Ref.  1  gives  the 
following  steady  state  estimates: 

AT(surface)  =  50°C 

AT( Interior ,  Y  =  1.6  cm)  =  70°C  (5) 

The  closeness  of  the  results  Indicates  that  the  basic  1-D  model  Is  sound. 

Now,  the  terms  W,  k,  and  h  will  be  examined  to  see  how  pad  temperatures  might 
be  reduced. 
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Here,  the  loading  of  the  pad  by  the  road  wheel  Is  modeled  as  simple 
one-dlmenslonal  compression,  and  the  pad  rubber  as  a  simple  linear 
visco-elastic  material  (small  strain).  The  heat  Input  per  unit  volume  per 
load  cycle  can  be  expressed  as*: 


(W/cps)  =  wtan(i)*(Elastlc  Strain  Energy) 

2 

2 

=  irtan(6)cr_,  =  1.1486  ps  1 ,  for  (6) 

4  E 

tan(6)  =  damping  factor  (=  .09  for  SBR) 

E  =  elastic  dynamic  modulus  (1040  psl  for  Interatlonal  Rubber  Hardness 
Degrees  (IHRD)  =  70) 

<j  =  average  compressive  stress  (=  130  psl  for  an  M60  tank) 
cps  =  loading  cycles  per  second 

The  average  number  of  compressive  load  cycles  for  a  given  time  can  be 
expressed  as 

V*t*N 

N  =  Number  of  load  cycles  =  RW 

c  L 

T  (7) 

V  =  velocity  of  vehicle  (20  mph  =  29.3  ft/s) 
t  =  time 

Nrw  =  number  of  road  wheels  (6) 
l.j.  =  length  of  track  (=  40  ft) 

then, 

N  V*N 

cps  =  __c_  =  "RW  =  4.4  (8) 

t  L 

W  =  (W/cps)*cps  =  5.05  ps1*s~l 

=  3.49*105  erg*cm~3»s_1 

This  Is  close  to  the  value  given  In  the  literature5. 
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Equation  6  shows  that  the  means  available  to  lower  W  are  to  lower  the  damping 
factor  and  to  raise  the  elastic  modulus.  The  mode  of  reinforcement  shown  In 
Figure  5-2  will  enforce  a  zero  strain  condition  In  the  wire  (X)  direction. 

For  simple  vertical  ( Y)  loading,  and  zero  lateral  (Z)  stress,  simple 
elasticity  gives  the  following: 

E  =  elastic  modulus 

v  =  Poisson  ratio  (=:  0.5  for  rubber) 

c  ,  c  =  normal  strains  In  the  X  and  Y  directions 
x  y 

o  ,  a  ,  a  =  normal  stresses  In  the  X,  Y  and  Z  directions 
x  y  z 


c  =0=l(d-vo),d  a  VO 
x  — '  x  y  x  y 

2 

c  =  1(<JV  -  VO  )  =  lo  (1  -  V  ) 

y  e  y  x  e  y 

°y  =  ^r'y-  fEey  (,0) 

1-v  J 

The  effective  modulus  In  the  Y  direction  Is  raised  by  1/3  due  to  the 
constraint  of  the  reinforcement.  Because  cx  =  oz  =  0,  all  of  the  work 
Is  still  done  by  oy.  Equation  6  (for  uniaxial  loading)  thus  still 
applies,  but  with  the  effective  modulus  as  shown  above.  Thus,  W  Is  reduced  by 
25%.  Reinforcement  In  a  second  direction  would  raise  the  pad  stiffness  to  the 
order  of  the  bulk  modulus  of  the  rubber  (~  100,000  psl),  which  would 
probably  Impair  the  suspension  functions  of  the  pad.  1-D  reinforcement  should 
be  acceptable,  since  the  stiffness  remains  of  the  same  order  as  when  no 
reinforcement  Is  present. 

The  damping  factor  can  be  changed  by  changing  the  material.  While  tan(5)  Is 
taken  as  0.09  for  SBR,  a  value  of  0.06  Is  usually  reported  for  natural  rubber 
(NR)  compounds.  Compounding  changes,  such  as  adding  different  amounts  or 
types  of  carbon  black,  will  directly  change  the  elastic  modulus. 

Metallic  wire  reinforcement  will  directly  affect  the  thermal  conductivity  In 
the  pad.  For  small  wire  volume  fractions,  the  effect  transverse  to  the  wire 
direction  can  be  Ignored.  In  the  wire  direction,  thermal  conductivity  at 
steady  state  will  be  a  simple  volume  average  of  the  values  for  the  wire  and 
the  rubber. 

A  specific  example  Is  as  follows:  The  wire  Is  carbon  steel,  0.032  Inch 

(0.8  mm)  diameter,  spaced  0.25  Inch  (6.4  mm)  apart  In  a  square  array,  oriented 

unlaxlaily. 

kwire  =  28  Btu*ft/(ft2*hr*°F)  =  4.85*10*  erg*cm/( cm2*s*°C) 

^rubber  =  3-2*10  erg*cm/(cm  *s*°C) 
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wire  volume  fraction  =  0.0129 

k  -  VW  *  <’  -  V*krubber 

k  =  9.4*10*  erg*cm/(cm  *s*”C)  (11) 


This  Is  a  significant  Increase  over  the  value  for  the  rubber.  It  cannot, 
however,  be  used  directly  In  the  1-D  formulation  (equation  3),  because  the 
wire  direction  Is  not  the  direction  of  the  minimum  pad  dimension  (the  Y 
direction).  Later,  2-D  numerical  results  will  be  presented  which  use  this 
Increased  thermal  conductivity  In  the  wire  direction. 

Surface  heat  transfer  can  be  Influenced  by  the  free  stream  temperature  Tf, 
and  by  the  heat  transfer  coefficient  h.  Tf  Is  effectively  the  temperature 
of  the  heat  sink  Into  which  energy  Is  transferred  from  the  body  In  question, 
h  Is  a  property  of  the  free  stream  material,  velocity,  manner  of  flow  (laminar 
vs  turbulent),  the  size  of  the  body,  and  other  factors.  For  Instance,  h  for  a 
free  stream  liquid  Is  usually  higher  than  for  a  free  stream  gas,  all  other 
factors  (1e,  free  stream  velocity)  being  equal.  The  form  presented  for  a  six 
Inch  body  length  and  100°F  free  stream  air  with  laminar  flow  over  the  pad 
Is7 


_  1  _5 

h  =  0.0157  Btu°F"Xhr  2  ft2  * 


h( 20  mph)  =  5.11  Btu/( ft2hr°F) 

=  2.90*10*  erg/(cm2s°C) 


(12) 


The  functional  form  shows  the  Influence  of  the  free  stream  velocity 
The  computed  value  at  20  mph  Is  somewhat  higher  than  that  used  In  the  1-D 
model.  For  the  2-D  numerical  studies,  a  compromise  value  Is  used  for  20  mph 
travel : 

h  =  2 . 35*1 04  erg/(cm2s°C)  (13) 

Equation  12  shows  that  Increasing  free  stream  velocity  can  Increase  surface 
heat  transfer;  this  suggests  that  a  blower  system  might  be  useful.  Another 
way  to  Increase  h  Is  to  Increase  the  effective  surface  area.  For  Instance, 
short  slots  on  the  X  surface  (normal  to  the  wire  direction)  could  double  the 
effective  area  on  that  surface,  thus  doubling  h  In  the  X  direction. 

5.2.3  Mechanical  Loads  From  Sharp  Asperities 

The  type  of  loading  which  results  from  running  over  sharp  asperities  ( 1 e , 
rocks,  debris)  can  be  seen  deduced  from  a  calculation8,  which  Is 
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sunmarlzed  In  Figures  5-4  and  5-5.  (The  calculation  Is  for  a  unit  thickness, 
plane  stress.)  The  angle  y  arises  from  the  action  of  the  road  wheel,  and 
causes  tensile  stresses  In  parts  of  the  pad.  The  orientation  of  this  tensile 
stress  helped  motivate  the  choice  of' the  direction  of  reinforcement. 

5.2.4  Stress  In  Reinforcement 

For  the  reinforcement  to  work  properly,  the  wires  must  not  break,  and  must 
maintain  adhesion  with  the  rubber.  Simple  estimates  of  the  tensile  stress  In 
the  wires  and  the  shear  stress  at  the  wire-rubber  Interface  are  given  In 
Figures  5-6  and  5-7.  A  high  value  of  compressive  stress  Is  used  to  model  the 
effect  of  very  uneven  terrain.  The  results  show  that  a  high  strength  wire  Is 
preferable.  Steel  tire  cord  Is  more  than  adequate,  while  aluminum,  copper,  or 
brass  wire  would  require  higher  volume  fractions. 

5.2.5  Fracture  Mechanics  Approach  To  Debonding 

In  the  above  discussion,  debonding  Is  approached  by  estimating  the  shear 
stress  at  the  bond  Interface.  A  better  approach  Is  to  use  a  fracture 
mechanics  formulation9*10.  The  basis  of  the  formulation  Is  a 
Griffith-type  energy  model  In  which  energy  for  debonding  comes  from  the  net 
loss  In  the  sum  of  the  potential  energy  of  the  load  and  the  strain  energy  In 
the  rubber.  (The  reinforcement  Is  assumed  rigid,  hence  storing  no  strain 
energy.)  Energy  conservation  can  be  written  as: 

^(Surface  energy  +  Potential  energy  +  Strain  energy)  =  0  (14) 

(4  Indicates  a  differential) 

One  development  assumes  a  constant  displacement  of  the  loading  device;  hence, 
the  potential  energy  does  not  change,  and  strain  energy  falls  as  the  debonded 
or  torn  area  grows  due  to  unloading  of  material  near  the  new  free 
surface9.  The  other  development  assumes  a  constant  load10.  In  that 
formulation,  potential  energy  drops  and  strain  energy  rises  as  debonding 
progresses . 

An  alternate  failure  mechanism  of  tearing  through  rubber  rather  than  debonding 
reinforcement  Is  also  explored9*10.  This  Is  an  unlikely  failure 
mechanism  for  reinforced  track  pads.  Analysis  will  concentrate  on  the 
debonding  problem. 

A  test  method  has  been  developed  for  measuring  the  adhesive  fracture  energy 
per  unit  area  (Ga)  needed  to  debond  a  wire-rubber  Interface10.  The  load 
Is  fixed  and  Ga  Is  Independent  of  crack  length,  so  when  debonding  starts, 

It  Is  total,  and  the  test  piece  (Figure  5-8)  falls  completely.  This  value  of 
Ga  represents  one  cycle  to  failure,  but  as  discussed  before,  track  pads  must 
withstand  many  cycles  to  failure.  Existing  data9  can  be  used  to  relate 
lower  values  of  Ga  to  debond  growth  per  cycle.  Hence,  some  estimate  of 
fatigue  life  can  be  made. 
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Figure  5-4.  Track  shoe  moving  over  surface  asperity 


cr  =#  ■  cos  (  y  +  e  > 

r  Trr 

crx  =  crrC0S2  (  0  )<  cry  =  crrSIN2  (  0  ),  =  CTrSIN  (  0  )  COS  (  0  ) 

ON  SURFACE,  AT  0  =  +  ,  <j  =  t2E  SIN  (  t  ) 
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(TENSILE  STRESS  AT  SURFACE) 
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Figure  5-5.  Mathematical  form  of  the  stresses 
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Figure  5-6.  Uniaxial,  fiber  loading  in  pad  under  compression 
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Figure  5-7.  Shear  stress  at  fiber-rubber  interface 
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Referring  to  the  wire  pullout  specimen  shown  In  Figure  5-8,  the  energy  balance 
(Equation  14)  can  be  written 


A[G  *2*ac 

a 


Fp/3  ♦  JqFcI/3]  =  0, 


0  =  C*AL 
L 


(15) 


0  Is  the  displacement  of  the  wire  which  occurs  when  the  newly  debonded  length 
c  stretches  by  the  strain  Al/l0.  The  first  term  In  brackets  Is  the  new 
surface  energy  created  when  the  debond  grows  by  length  c.  The  second  term  Is 
the  drop  In  potential  energy  of  the  applied  load  due  to  the  dlsplacment  0. 

The  third  term  Is  the  gain  In  strain  energy  due  to  stretching  of  the  rubber 
adjacent  to  the  debonded  length  c  by  the  strain  AL/L0.  Applying  the 
differential  to  equation  15  gives 


G,  =  1  [F  AL  -  JnLFd(AL)] 

2iraL  P  U 

0 


(16) 


The  term  In  brackets  Is  a  nonlinear  complementary  energy  function,  as  shown  In 
Figure  5-9.  For  linear  stress-strain  behavior.  Equation  16  can  be  simplified 
as  follows: 


f-  a 

A  i  A  -  ira  ,  a  =  p ,  c  *  4k  = 

P  A  L  E 

0 

jJLFd(AL)  =  FpAL  (17) 

2 


2 


4irAaE 


This  Is  the  form  derived  In  the  literature10. 

As  discussed  before,  Ga  from  Equation  16  or  17  Is  a  surface  energy  term  for 
failure  In  one  loading  cycle.  When  the  applied  load  Is  less  than  the  failure 
load  (Fp),  the  test  piece  can  withstand  more  than  one  cycle  to  failure.  The 
following  example  shows  how  the  number  of  cycles  to  failure  for  a  less  than 
critical  load  could  be  estimated. 

Consider  a  wire  pullout  specimen  with  linear  elastic  behavior  (Equation  17) 
and  with  the  following  conditions  at  failure: 
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F  =  32  lbf 
P  a 

A  =  0.25  In 

L  =  1  Inch,  Al  =  0.7  Inch 
E  =  183  psl 
a  =  0.016  Inch 

G  =  111.3  lb/in  =  19.5  kJ/m  (Equation  17) 
d 

This  Is  typical  of  the  values  of  Ga  reported  In  Reference  10  for  brass 
coated  steel  wire  In  rubber.  It  represents  one  cycle  to  failure. 

Now, consider  the  case  of  cyclic  loading  with  a  maximum  force  of  8  lb.  Using 
Equation  17,  this  gives  an  available  debonding  energy  per  unit  area  G  =  1.2 
kJ/m2.  The  data  from  Reference  9  Is  presented  In  Figures  5-10  and  5-11, 
which  give  the  debond  growth  per  cycle  as  a  function  of  the  available 
debonding  energy.  The  data  for  reinforced  rubber  Is  slightly  to  the  left  of 
that  for  rubber  alone,  Indicating  that  the  bond  energy  Is  somewhat  less  than 
the  tear  energy  In  the  rubber  Itself.  Using  the  dashed  line  In  Figure  5-11, 
the  debond  growth  per  cycle  Is  roughly  40  nm  (nano-meters),  or  4  x  10~* 
cm.  Thus,  for  the  debond  to  grow  by  1  cm,  250,000  load  cycles  would  be 
required. 

A  calculation  similar  to  the  one  above  can  be  carried  out  for  the  reinforced 
track  pad  configuration,  shown  In  Figure  5-12  with  the  pad  Idealized  as  a 
rectangular  block.  Since  the  loading  Is  through  a  steel  plate,  both  bonded 
and  debonded  sections  of  the  pad  see  the  same  strain.  The  rubber  Is  Idealized 
as  linear  elastic,  and  a  small  strain  elastic  modulus  Is  used.  As  derived  In 
Equation  10,  the  effective  modulus  In  the  bonded  region  Is  taken  as  4/3  the 
modulus  of  the  rubber  alone.  Bonded  and  debonded  regions  thus  see  the  same 
vertical  strain,  but  different  stresses.  Recognizing  that  the  debonded  length 
c  Is  the  variable.  Equation  14  can  be  written  as 


_d(Surface  Energy)  =  ^(Potential  Energy  +  Strain  Energy)  (19) 
dc  dc 

For  linear  elastic  behavior,  the  gain  In  strain  energy  equals  half  of  the  loss 
In  potential  energy,  as  Is  the  case  In  Equation  17.  The  terms  In  Equation  19 
are  computed  as  follows: 


F  -  (Lx-c>LziER‘,  *  cHzVy 


■  LZER<iLx 

o 


c) 

3 


Potential  Energy  (P.E.)  = 

p,.  .  -E\ 


-F*4Y  = 


_FcyLy 


Lz¥JLx  -  £>’ 


(20) 
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Figure  10.  Comparison  of  averaged  crack  growth  rates  in  laminates  of  standard  construction  (o)  with  resuf 
from  independent  measurements  on  tensile  strip  test-pieces  of  the  ply  rubber  (•).  The  bars  indicate  th 
range  of  the  results  for  the  laminates. 


Figure  5-10.  Comparison  of  Averaged  Crack  Growth  Rates  (Ref.  9) 
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Averaged  crack  growth  rates  lor  laminates  ol  various  constructions.  The  line  is  the  same  ns  that 
draw  n  through  the  independently  measured  nark  growth  characteristics  shown  in  li^ure  ft*. 


Figure  5-11.  Debond  Growth  Per  Load  Cycle  (Ref.  9) 
Scales  as  in  Figure  5-10,  Dashed  Line 
Added  Here. 
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Surface  Energy  (Su.E.)  =  GwDcN^ 

Nw  =  number  of  wires  In  the  block,  0  =  wire  diameter 

2 

D  N 

wire  volume  fraction  =  V,  =  *  _ w 

r  4  L  L 
Y  z 

Su.E.  =  GtrDc*4VfLyLz 

2 

*0 

As  discussed  before, 

d_(Stra1n  Energy)  =  -I  d  (P.E.) 
dc  2  dc 

Equation  19  then  becomes 


G«0*4VflyLz  .  _ T\l 


irD 


6LzER<fx 


c)‘ 

3 


G  = 


OF 


24-VfLzER<!Lx  -  |)2 


Typical  numbers  are 


Wire  diameter  D  =  0.032  Inch 

Load  F  =  5000  lbf 

Wire  Volume  fraction  V^  =  0.013 

L  =8  Inches 
z 

rubber  elastic  modulus  E  =  1000  psl  (small  strain) 

K 

L  =4.5  Inches 
x 

debonded  length  c  =  1  Inch 

G  =  7.07  lb/in  =  1.24  kJ/m* 


(21) 


-  (22) 


(23) 


Using  Figure  5—11  (dashed  line),  this  corresponds  to  6*10"*  cm  of  debond 
length  per  load  cycle. 

This  can  be  translated  Into  a  relation  between  mileage  and  debonded  length  as 
follows: 
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track  length  L  =  40  ft 
number  of  road  wheels  =  6  , 
distance  travelled  X  =  1000  miles 

number  of  load  cycles  N  =  *  ^RW  =  7.92*10 

4 

c  =  6*1 0-6  cm/cycle*Nc ,  =  4.75  cm,  =  1.87  Inches  (24) 


This  Indicates  that  some  significant  debonding  will  take  place  over  that 
distance.  Equation  22  Indicates  how  Improvements  could  be  made.  Wire 
diameter  could  be  made  smaller,  and  wire  volume  fraction  could  be  Increased. 
Going  to  a  finer  wire  size  (.016  Inch  diameter)  or  a  higher  Vf  (.026)  would 
cut  G  by  a  factor  of  2  (to  0.62  kJ/m2) .  The  debond  length  per  cycle 
(dc/dNc)  from  Figure  5-11  would  then  be  6  x  10~?  cm;  one  tenth  of  the 
previous  value.  Debonded  length  would  then  be  0.187  Inch. 

Several  approximations  were  made  so  that  this  analysis  could  proceed;  It  Is 
useful  to  discuss  their  Impact  on  the  results.  First,  the  effect  of  friction 
on  the  debonding  process  was  not  considered.  Since  the  pad,  and  the  bonds, 
are  loaded  In  compression,  this  might  significantly  decrease  dc/dNc.  Second, 
temperature  effects  may  Increase  dc/dNc  as  the  pad  temperature  Increases. 

(See  section  on  material  development,  and  Appendix  B) .  Third,  the  use  of  the 
data  from  Figure  5-11  In  this  manner  Involves  approximations  due  to  both  the 
scatter  In  that  data  and  the  extrapolation  to  a  different  test  geometry.  The 
results  still  retain  validity,  but  the  trends  may  be  more  accurate  than  the 
absolute  numerical  results.  For  Instance,  the  example  In  the  previous 
paragraph  shows  that  fatigue  life  Is  very  sensitive  to  the  bonded  area,  as 
reflected  In  the  wire  size  and  wire  volume  fraction.  If  a  given  choice  of 
wire  diameter  and  volume  fraction  gives  marginal  results  In  testing,  modest 
changes  In  these  parameters  could  result  In  large  Improvements  In  performance. 


5.2.6  2-D  Time  Dependent  Thermal  Solution  Method  via  Finite  Differences 


The  track  pad  size  Is  roughly  2  Inches  x  4.5  Inches  x  10  Inches.  The  more 
Important  temperature  variations  will  be  seen  In  the  short  dimensions,  so,  for 
the  2-D  model,  the  long  dimension  Is  omitted.  The  grid  model  Is  shown  In 
Figure  5-13.  X  and  Y  were  used  as  the  horizontal  and  vertical  directions.  A 
6x6  grid  was  found  sufficient  to  accurately  model  typical  test  cases,  as 
were  found  In  References  4  and  11.  (These  were  1-D  and  2-D  problems, 
typically  solved  by  expansion  techniques). 

The  governing  equation,  analagous  to  the  first  of  Equations  2,  but  with  time 
dependence  and  two  spaclal  dimensions  retained,  Is 


pcaT  =  k 
at 


2  2 

a  T  +  k  a  T  + 

*2  y  2 

ax  3Y 


w 


k^  =  k  In  X  direction,  Increased  by  presence  of  wire 
k  =  k  In  Y  direction,  =k  of  rubber 

y 


(25) 
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44  INCHES 


(1,1) 


SYMMETRY 


Two  types  of  boundary  conditions  will  be  considered.  The  first,  a  convective 
boundary  condition,  Is  analagous  to  the  second  of  equations  2: 


¥  -  ¥  (rs  '  V 

3n  k 

n 


(26) 


n  =  outward  normal  direction. 

This  would  typically  apply  to  a  body  sitting  In  an  air  stream.  The  second 
type  of  boundary  condition  Is  a  specified  surface  temperature,  and  would  apply 
to  a  body  held  In  a  heated  mold: 

TS  =  Tm  (27) 

The  finite  difference  time  step  requires  a  finite  difference  form  of 
Equation  25.  Where  (I,J)  represent  node  points  (I  varies  In  the  X  direction, 
J  In  the  Y  direction)  In  Figure  5-13,  and  AX,  AY  are  the  distances  between 
node  points,  the  typical  finite  difference  form  of  the  second  derivative  Is 

a2T  =  T( 1+1 , J)  -  2T( I , J)  +  T( 1-1 , J) 

2  2 

ax  (ax) 


This  Is  modified  at  the  line  of  symmetry  (I  =  1)  and  at  the  boundary  ( I  = 

NI).  At  the  line  of  symmetry,  1=1,  and  T(I-1,J)  will  be  equal  (by  symmetry) 
to  T(I+1,J).  Therefore 

all  =  2 [ T ( 2 , J )  -  T( 1 , J ) ] 

ax2  (AX)2  (29) 

1=1 


At  the  boundary,  first  consider  the  condition  given  by  Equation  26.  The  first 
derivative  at  the  surface  can  be  put  Into  difference  form: 


il  =  T(NI  +  1  .  J)  -  T(NI-1.J).  =  [T(NI ,  J)  -  Tf] 

aX  2 ( AX)  k  (30) 

I  =N  I  X 


T(NI+1,J)  Is  a  fictitious  value,  lying  outside  of  the  body.  It  Is  useful 
because  It  can  be  computed  from  Equation  30,  and  then  used  to  compute  the 
second  derivative: 


T(NI  +  1 , J)  =  TQut  =  T(NI-1 , J)  -  2( AX)_x[T(NI , J)  -  T  ] 

k  r 

x 


(31) 
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A  similar  development  applies  for  a2T/aY2. 

When  Integrating  forward  In  time,  the  time  step  must  be  chosen  small  enough  to 
maintain  accuracy.  Some  guidance  can  be  found  In  Reference  11.  Consider  the 
1-D  difference  form  of  Equation  25: 

PCAI  =  krT(I+l)  -  2T(  I )  +  T(I-D] 

6t  (M)’  («> 


From  dimensional  considerations,  one  can  write 


At  Is  order  of  [(AX)  gc] 

k 

2  _ 

=  constant*! AX)  ,  0  =  k_ 
FT  Pc 


(33) 


The  constant  Is  sometimes  taken  as  0.5  for  a  special  solution  method11. 

Here,  It  Is  taken  Initially  as  0.2,  At  Is  calculated  for  the  X  and  Y 
directions,  the  smaller  value  chosen,  and  further  reduced  to  some  convenient 
Integral  value.  (In  calculations  presented  later,  At  Is  25  seconds,  reduced 
from  the  33.39  seconds  calculated  for  the  Y  direction). 

5.2.7  Curing  Model 

So  far,  the  formulation  has  concentrated  on  modeling  the  hysteresis  heating  of 
fully  cured  rubber  track  pads.  Because  the  formulation  calculates  temperature 
throughout  the  pad  as  a  function  of  time,  and  because  a  fixed  temperature 
boundary  condition  can  be  used,  the  curing  process  Itself  can  be  modeled. 

This  Is  Important  because  thick  rubber  parts  are  difficult  to  cure  evenly,  due 
to  the  low  thermal  dlffuslvlty  of  rubber,  but  good  performance  requires  a 
fairly  even  cure.  The  difficulty  arises  because  the  outside  of  the  block, 
being  close  to  the  mold  walls,  usually  sees  a  more  severe  temperature-time 
history  than  the  Inside  of  the  block. 

A  common  approximation  used  to  model  the  dependence  of  rate  of  cure  on 
temperature  Is  the  temperature  coefficient  of  vulcanization  (TC),  and  Is  the 
factor  by  which  the  rate  of  cure  Increases  for  a  10°C  Increase  In  temperature 
over  some  reference  temperature  (T^p)11.  A  standard  value  from  the 
literature  Is  TC  =  2.  The  equation  used  here  to  describe  this  process  Is 


d(  SC) 
dt 


(TC) 


[0.1 (T-Tref)  ] 


(34) 
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where  SC  Is  state  of  cure.  As  a  simple  example,  consider  an  Ideal  cure  which 
Is  27.5  minutes  at  290°F  (143. 3°C).  Taking  143. 3°C  as  the  reference 
temperature.  Equation  34  gives  rates  of  curing  for  differing  temperatures. 

The  actual  state  of  cure.  In  equivalent  time  at  T  =  T^p,  can  be  expressed  as 


SC  =  I*d(SC)dt  =  J*(TC) 
0  dt  0 


[0.1(T(t)-TREF)] 


This  Is  easily  added  to  the  finite  difference  formulation,  which  updates  the 
temperature  at  each  node  as  time  progresses.  A  different  curing  relation, 
such  as  a  different  equation  or  tabular  data,  could  also  be  used. 


5.3  Material  Development 


5.3.1  Material  Selections 


A  range  of  rubber  compounds  was  considered,  based  on  physical  properties  and 
adhesive  bonding  to  reinforcement.  A  series  of  test  compounds  were  devised 
based  on  Bergstrom's  baseline  SBR12,  and  Included  both  elemental  sulfur 
(conventional)  and  EV  (sulfur  donor,  efficient  vulcanization)  cure  systems. 
Compounds  with  conventional  cures  can  bond  to  bare  brass  or  bronze  coatings, 
with  no  adhesive  needed.  EV  cure  systems  sometimes  give  better  physical 
properties,  and  can  also  be  less  sensitive  to  overcuring.  (The  Inslde-of  the 
pad  can  be  more  fully  cured  without  causing  reversion  (property  degradation) 
near  the  surface).  However,  EV  cure  systems  do  not  give  bonding  to  brass  or 
bronze  coatings,  so  an  adhesive  must  be  used.  If  the  decision  Is  made  to  use 
an  adhesive,  then  brass  or  bronze  coatings  are  no  longer  a  constraint  with 
regard  to  bonding.  (They  are  still  useful  for  corrosion  protection,  as  was 
found  when  test  specimens  were  prepared.)  Other  rubber  compounds  tested  were 
an  all  NR  compound  with  conventional  cure,  a  compound  from  a  companion 
program,  and  a  TACOM  standard  triblend  (14A,  MIL-T-1 1 891 C( AT) )  conventional 
cure  compound. 


Some  rubber  compounds  used  In  this  study  are  given  In  Tables  5-1  and  5-2. 

Table  5-1  shows  compounds  made  expressly  for  this  program  (Avon  Custom  Mixing, 
Holbrook^  Mass).  Table  5-2  shows  a  compound  made  for  the  companion 
program13,  which  was  tested  In  this  program.  The  14A  compound  Is  detailed 
In  the  MIL  spec.  Samples  were  supplied  by  Avon  Custom  Mixing,  and  Goodyear 
(custom  compounded  by  Avon,  from  stock  by  Goodyear).  A. .conventional  cure 
compound  of  unknown  composition  was  also  supplied  by  Avon  for 
experimentation.  Some  Important  features  of  the  compounds  are  summarized 
below. 

Designation  Rubber  Type  Cure  Type 


20A 

SBR 

conv.(h1gh  sulfur) 

20B 

SBR 

EV  (low  sulfur) 

20C 

SBR 

conv.  (high  sulfur) 
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200 

SBR 

conv.  (med. 

sulfur) 

20E 

NR 

conv.  (high 

sulfur) 

A54 

'  SBR/BR/NR 

conv.  (med. 

sulfur) 

14A 

SBR/BR/NR 

conv.  (med. 

sulfur) 

Here,  high  sulfur  means  3.5  phr,  which  Is  a  standard  recommendation  for 
adhesion  to  brass. 

Four  types  of  wire  were  considered:  solid  brass  wire,  brass  coated  steel  wire 
(tire  wire)  bronze  coated  steel  wire  (hose  wire),  and  a  steel  wire  with  a 
thick  copper  coating.  The  latter,  whose  trade  name  Is  CopperplyR,  has  30% 
of  the  thermal  conductivity  of  solid  copper,  yet  retains  the  strength  of 
steel,  since  It  Is  roughly  80%  steel  by  volume.  (Samples  of  tire  wire,  hose 
wire,  and  CopperplyR  were  obtained  from  National-Standard  Company,  Niles 
Mich.). 

Samples  of  several  types  of  metal-rubber  adhesives  were  obtained  from  Lord 
Chemical  Corporation,  Erie,  Pa.  These  can  be  applied  to  the  wire,  which  Is 
then  Inserted  Into  the  uncured  rubber.  The  adhesive  then  co-cures  with  the 
rubber.  Before  curing,  the  adhesive  Is  not  bonded  to  the  wire,  so  the  surface 
coat  Is  fragile.  The  adhesive  of  choice  was  ChemlokR  252,  a  solvent  based 
adhesive.  It  was  chosen  because  It  appeared  to  be  the  easiest  to  apply,  and 
gave  a  relatively  durable,  flexible  coat  In  the  uncured  state.  (Water  based 
adhesive  coats  were  brittle  when  dry).  To  clean  a  copper  surface  In  „ 
preparation  for  this  adhesive,  the  following  was  used,  based  on  the 
manufacturer’s  recommendation;  a  one  minute  etch  In  a  25%  by  weight  solution 
of  ammonium  persulfate,  followed  by  a  water  rinse.  Initial  cleaning  was  an 
acetone  wash  and  an  air  or  towel  dry.  Note  that  an  etch  was  not  required  for 
the  wires  with  brass  or  bronze  surfaces  because  these  were  resistant  to 
oxidation.  The  adhesive  was  applied  by  dipping,  then  dried  In  an  oven  at  low 
heat.  When  necessary,  the  adhesive  was  diluted  with  solvent  to  give  a  thin 
coat,  as  per  Instructions  from  the  supplier.  Coated  wires  were  sometimes 
stored  for  several  days  before  the  specimens  were  molded  and  cured. 

5.3.2  Adhesion  Tests 

A  series  of  adhesion  tests  were  conducted  as  described  In  Equations  14 
through  18.  The  test  specimens  were  rubber  blocks,  0.5  Inch  x  0.5  Inch  x  3 
Inches  long,  with  wires  Inserted  1  Inch  Into  the  block,  along  the  centerline, 
from  each  end  (see  Figure  5-8).  Samples  were  prepared  by  laying  up  the 
uncured  rubber  and  the  wires  In  a  special  mold,  and  curing  In  a  heated  press. 
Typical  cures  were  30  to  60  minutes  at  290°F  In  the  press,  then  cooling  under 
tap  water  and  room  temperature.  Samples  were  tested  using  an  Instron  testing 
machine,  with  a  load  rate  of  50  pounds  per  minute.  As  discussed  before,  the 
proper  loading  method  for  the  test  Is  prescribed  force  (dead  load),  rather 
than  prescribed  displacement. 

The  adhesion  test  results  are  summarized  In  Table  5-3.  Ga  Is  calculated 
using  the  nonlinear  formulation  (Equation  17),  and  a  linear  small  strain 
approximation  using  Equation  18  with  an  engineering  stress-strain  secant 
modulus.  This  gives  the  following: 
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Table  5-1  Rubber  Formulations  (PHR) 


SBR  1500 

100 

100 

N300  HAF 

45 

45 

ZnO 

5 

5 

Stearic  Acid 

2 

2 

UOP  88 

3 

3 

Neozone  D 

1 

1 

Hellozone  Wax 

1 

1 

Process  Oil 

3 

3 

Sulfur 

3.5 

0.2' 

Amax 

0.8 

- 

Sulfsan  R 

0.5 

- 

Methyl  Tuads 

- 

0.5 

Morphax 

- 

2.5 

Santocure  NS 

- 

- 

NR 

Zinc  Stearate 


Santocure  MOR 


[i 


Table  5-2  Long-Life  Tank  Tread  (REF.  13)  Formulation  A54 

A54 


Phil Iprene  1609 

40.6 

Cls-4  13501  0-  Cls-4  1352 

25.8 

Pale  Crepe 

60 

Saf  Black 

42 

Zinc  Oxide 

3 

Stearic  Acid 

2 

Sulfur 

2 

Santocure 

1.5 

Thermoflex  A 

1 

Plccopale  100 

3.5 

Santoflex  AW 

1.5 

U0P  88 

5 

Hellozone  Wax 

1 

1C1s-4  Polybutadiene  master  batch  (100  parts  C 1  s -4  1203,  80  parts 
Ph 1 1  black  I  and  35  parts  Phi  1  rich  5). 

2C1s-4  Polybutadiene  master  batch  (100  parts  C 1  s -4  1203,  90  parts 
Phllblack  I  and  50  parts  Phi  1  rich  5). 


-34- 


Table  5-3:  Wire  Adhesion  Tests  (all  wire  diameters  0.032  Inch  except  as  noted) 
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Table  5-3:  Wire  Adhesion  Tests  (Cont'd) 
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E  =  a  =  FPL0. 
c  AAL 


G( 1  Inear)  =  FpAL 
4*AL 


The  force-displacement  test  records  are  shown  In  Figures  5-14  through  5-22. 
They  Include  relatively  poor  results  and  very  strong  results.  The  curves  show 
a  distinct  nonlinearity,  especially  at  the  higher  displacements  (At).  In 
all  cases,  L0  =  1  Inch,  so  the  engineering  strain  equals  Al  In  Inches. 

Several  conclusions  were  drawn  from  the  test  results,  some  of  which  were 
supported  by  consultation  with  suppliers.  These  were: 

1)  Brass  coated  steel  wire  can  give  a  higher  adhesion  than  solid  brass  wire 
(Tests  15,  16  and  21  vs.  test  20).  The  supplier  of  the  wire  attributed 
this  to  diffusion  In  the  solid  brass  wire,  which  changes  the  surface  com¬ 
position.  This  cannot  occur  In  the  coated  wire,  since  the  brass 
components  cannot  diffuse  Into  the  steel. 

2)  The  adhesive  worked  well  for  the  EV  cure  rubber,  Increasing  Ga  for  a 
given  rubber-wire  combination  by  a  large  margin  (tests  1-5  vs.  test  13). 

A  surface  etch  Is  not  needed  with  brass,  and  may  even  reduce  adhesion,  but 
Vs  needed  to  clean  the  oxide  layer  from  copper  (tests  13,  14,  18,  19). 

3)  The  best  adhesion  was  obtained  for  a  high  sulfur,  conventional  cure, 
natural  rubber  compound,  with  brass  coated  steel  wire  (test  20). 

4)  For  the  TACOM  standard  (compound  14A),  adhesion  to  the  bronze  coated  wire 
was  roughly  twice  that  obtained  for  the  brass  coated  wire  (tests  22-24). 
The  supplier  commented  that  this  difference  In  adhesion  Is  a  function  of 
the  rubber  compound,  and  Is  not  unusual.  Better  results  might  have  been 
obtained  with  an  adhesive.  If  this  compound  were  to  be  used  to  make  steel 
wire  reinforced  pads,  a  bronze  coat  with  or  without  adhesive,  or  a  brass 
or  copper  coat  with  adhesive,  would  be  recommended  for  trial. 

5)  The  wrong  choices  can  lead  to  very  low  adhesion,  so  design  choices  should 
be  tested. 

According  to  the  supplier,  brass  coating  withstands  storage  better,  because  It 
Is  less  sensitive  to  moisture.  Bronze  coated  steel  wire  must  therefore  be 
stored  carefully.  In  a  low  moisture  environment,  before  use  as  reinforcement 
In  rubber. 

Some  descriptions  of  other  adhesion  tests  refer  to  a  tearing  of  rubber  rather 
than  a  debonding  of  the  Interface,  Indicating  that  the  bond  was  stronger  than 
the  rubber  Itself.  In  all  of  our  tests,  the  wires  pulled  out  with  clean 
surfaces,  Indicating  a  failure  of  the  bond  rather  than  the  rubber. 
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Figure  5-14  Wire  pullout  adhesion  tests 
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Figure  5-19.  Wire  puiiout  adhesion  tests 
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Equation  17  shows  that  for  equal  values  of  Ga,  the  rubber  with  the  higher 
modulus  will  have  the  higher  pullout  force.  This  could  be  an  Important 
consideration  when  comparing  adhesive  performance,  when  specimens  have 
different  geometries,  or  different  rubber  properties.  In  other  words,  Ga  Is 
a  property  of  the  adhesive  Interface,  and  Is  Independent  of  rubber  stiffness 
or  specimen  geometry. 

A  brief  comment  Is  Indicated  on  how  to  handle  different  measured  values  of 
Ga  In  terms  of  computing  number  of  cycles  to  failure  at  an  available  surface 
energy  G.  One  possible  method  Is  to  scale  G  by  the  ratio  of  some  reference 
value  of  Ga  to  the  actual  test  value.  A  typical  value10  of  Ga  Is  20 
kJ/m2,  so  set  Ga(ref)  =  20  kJ/m2.  Then, 


*  Ga(ref.) 

G  "  G*G  (meas.)  (37) 

a 

A  higher  measured  value  of  Ga  results  In  a  lower  value  of  G*,  which  then 
gives  more  cycles  to  failure. 

5.3.3  Effect  of  Temperature  on  Tire  Cord  Adhesion 

A  recent  telephone  Inquiry  to  the  Malaysian  Rubber  Bureau  on  this  subject 
resulted  In  the  following  unpublished  data  for  adhesion  of  brass  plated  steel 
tire  cord  to  rubber. 14 

Temperature  °C  Steel  Cord  Adhesion  (Newtons) 
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1 ,207 

50 
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60 
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70 
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80 

1  ,015 

100 
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120 

891 

140 

793 

The  NR  compound  used  In  this  study  consisted  of:  NR  100,  ZnO  5,  Zinc  ethyl 
hexanoate  2,  Flectol  H  0.5,  HAF  black  40,  oil  4,  Silica  10,  Resorcinol  2, 
Hexamlne  2,  Santocure  M0R  0.7  and  Sulfur  3.  The  cure  condition  used  for 
molding  the  steel  tire  cords  Into  the  rubber  Is  27  minutes  at  140°C.  The  test 
piece  wasjlmllar  to  that  used  In  ASTM  D2229,  with  12.5  mm  cord 
embedment15.  The  data  shows  that  there  Is  some  loss  of  adhesion  with 
temperature.  Temperature  also  affects  rubber  stiffness,  with  stiffness 
roughly  proportional  to  absolute  temperature  for  temperatures  between  -20°C 
and  +70°C.  Using  this,  plus  Equation  17,  the  change  In  G*  can  be 
written 
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(38) 


Ga^  =  T  =  absolute  temperature 

Ga(1)  [Pp(l)]2T(2) 


If  the  linear  dependence  of  stiffness  on  absolute  temperature  holds  to  140°C, 
the  reduction  In  Ga  from  the  value  at  room  temperature  Is  roughly  70%.  At 
70°C,  the  reduction  Is  roughly  40%.  This  level  of  change  Is  significant,  and 
shows  that  room  temperature  data  may  not  be  adequate  In  some  applications.  On 
the  helpful  side,  stiffening  due  to  Increased  temperature  will  decrease  the 
energy  available  for  debonding.  In  a  load  controlled  configuration. 

5 . 4  Summary 

Design  Ideas  are  presented  for  reinforced  tank  track  pads.  The  Ideas  are 
based  on  observations  of  TACOM  test  results  for  loading  of  a  track  pad  by  an 
asperity,  the  desirability  of  maintaining  rubber-like  compliance  In  a 
reinforced  pad,  the  advantage  of  reducing  pad  heat  build-up  from  cyclic 
loading,  and  the  need  to  maintain  a  strong  bond  between  the  reinforcement  and 
the  rubber.  An  outgrowth  of  the  analysis  were  some  observations  on  methods  of 
achieving  greater  uniformity  In  curing. 

The  main  features  of  the  suggested  baseline  design  are  steel-based  wire,  with 
brass,  bronze,  or  copper  coating,  aligned  In  the  direction  of  travel,  with  a 
volume  fraction  of  0.013  and  a  diameter  of  0.032  Inch.  See  Figure  5-23. 
Methods  are  presented  for  estimating  the  effects  on  performance  of  changes  In 
wire  volume  fraction  and  diameter. 

The  program  Included  thermal  and  mechanical  analyses,  plus  a  series  of 
adhesion  tests.  The  thermal  analysis  Included  a  simple  1-D  closed  form  model, 
and  a  numerical,  2-D,  time  dependent  model.  The  2-D  model  was  used  to 
simulate  both  heating  due  to  cyclic  loading,  and  the  curing  process.  Part  of 
the  mechanical  analysis  was  a  fracture  mechanics  based  treatment  of  debonding 
of  wire  reinforcement  as  a  function  of  load  magnitude,  configuration,  and 
number  of  loading  cycles.  A  test  method  was  chosen  which  yielded  adhesive 
surface  energy  per  unit  area.  In  the  testing  program,  good  results  were 
achieved  with  both  conventional  and  EV  cure  rubber  compounds.  Bonding  to  EV 
cure  rubber  required  an  adhesive.  It  was  very  sensitive  to  proper  surface 
preparation  and  adhesive  application.  Conventional  cure  compounds  could  bond 
well  to  brass  and  bronze  coated  steel  wire,  but  bond  strength  was  a  strong 
function  of  rubber  type  and  sulfur  content.  The  best  performance  was  measured 
for  a  conventional  cure,  natural  rubber  compound  with  high  sulfur  content. 
Tests  were  run  for  a  TACOM  baseline  triblend  (compound  14A),  with  brass  and 
bronze  coated  steel  wire.  The  bronze  coat  gave  the  better  performance. 
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APPENDIX  A  ..  FINITE  DIFFERENCE  CALCULATIONS 


A  computer  program  was  written  to  carry  out  calculations  using  the  analysis 
presented  In  the  previous  section.  Cases  were  run  for  heating  by  cyclic 
loading,  and  for  various  cure  cycles.  A  list  of  the  material  properties,  and 
graphs  of  temperature  and  cure  state  vs  time  for  cases  of  cyclic  loading  and 
curing,  are  contained  In  this  appendix.  For  the  cyclic  loading  cases,  steady 
state  values  of  temperature  were  calculated  using  a  series  solution  from 
Reference  4  for  a  rectangular  block  with  convective  boundary  conditions. 

Figures  A1  and  A2  show  temperature  at  the  pad  center  as  a  function  of  time  for 
various  loading,  reinforcement,  and  surface  conditions.  The  heating  rate  Is 
reduced  In  Figure  A2  to  account  for  the  stiffening  effect  of  the  wires.  As 
stated  before,  the  copper  coated  steel  wire  Is  a  commercially  available 
product,  which  Is  roughly  20%  copper  by  volume.  (Because  of  the  added 
complication  of  etching,  rinsing,  and  applying  an  adhesive  to  promote  bonding 
to  the  rubber,  this  product  was  not  as  attractive  as  the  brass  coated  steel 
wire).  In  all  of  the  calculations,  the  Initial  temperature  was  20°C.  When 
wires  were  present,  the  wire  volume  fraction  was  0.013.  The  relative  effect 
of  each  aspect  of  the  reinforcement  on  the  center  temperature  can  be  seen  by 
tabulating  the  steady  state  values: 


Baseline  T  =  85°C 

Stiffening  Alone  68 

r 

Stiffening  +  enhanced  kx  (steel)  61 

Stiffening  +  enhanced  kx(steel)  +  double  h^  56 

Stiffening  +  enhanced  k^  (copper  coated)  +  double  h^  51 


The  largest  single  effect  Is  the  stiffening,  since  this  directly  reduces  the 
rate  of  energy  Input  by  25%.  The  effect  of  the  wires  on  thermal  conductivity 
plus  the  effect  of  a  slotted  surface  on  surface  heat  transfer  are  somewhat 
smaller,  but  of  the  same  order.  A  reasonable  case  could  be  made  for  softening 
the  rubber  compound  so  that  a  reinforced  pad  would  have  the  same  compliance  as 
the  original  unreinforced  pad.  Then,  the  enhancement  of  kx  and  hx  would 
dominate.  Temperature  effects  alone  could  cut  the  tear  growth  rate  in  half 
(see  Figure  5-3).  Note  that  the  effect  on  kx  of  the  copper  coated  wire 
could  be  achieved  with  the  steel  wire  by  roughly  doubling  the  volume  fraction. 

The  next  set  of  graphs  represent  temperature  and  cure  state  vs  time  for 
several  cure  cycles  Imposed  on  steel  wire  reinforced  pads  (X  direction,  V f  = 
0.013).  The  target  value  of  cure  was  In  most  cases  27.5  minutes  at  290°F 
(143. 3°C),  which  Is  typical  of  what  Is  used  for  track  pad  rubber  compounds. 
Both  fixed  temperature  and  convection  boundary  conditions  are  used;  the  former 
represents  conditions  In  the  mold,  while  the  latter  represents  conditions 
outside  the  mold. 
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Figures  A3  and  A4  represent  a  simple  two  stage  cure,  In  which  a  block 
Initially  at  room  temperature  (20°C)  is  held  In  a  mold  at  270°F  (132. 2°C)  for 
65  minutes,  then  taken  out  to  cool  at  room  temperature.  The  temperature  plot 
(A3)  shows  that  the  center  of  the  pad  Is  at  deficit  In  the  mold,  but  retains 
temperature  longer  outside  the  mold.  The  cure  state  plot  (A4)  shows  that 
while  the  outside  of  the  pad  Is  reasonably  close  to  the  target  value,  the 
center  In  this  case  Is  very  much  undercured. 

Figures  A5  through  A8  show  two  different  three  stage  cure  cycles.  In  which  the 
block  Is  heated  In  an  oven  (convection  b.c.,  100°C)  before  being  put  Into  the 
mold  (temperature  b.c.,  132. 2°C).  When  removed  from  the  mold,  the  block  Is 
cooled  at  room  conditions  (convection  b.c.,  20*0) .  The  first  three  stage  cycle 

proved  to  have  too  little  time  In  the  mold;  not  even  the  outside  was  fully 

cured.  Time  In  the  mold  was  Increased  In  the  second  three  stage  cycle, 
resulting  In  a  cure  which  was  fairly  close  to  optimum  both  at  the  center  and 
the  surface  of  the  pad. 

Figures  A9  and  A10  represent  a  four  stage  cure  cycle,  in  which  an  oven 
post-heat  at  100°C  Is  added.  The  state  of  cure  Is  slightly  more  even  than 
that  of  the  previous,  three  stage  cure  cycle.  Time  In  the  mold  was  45  minutes, 
vs.  65  minutes  for  the  three  stage  cycle  (a  significant  reduction). 

In  the  previous  cases,  the  mold  temperature  was  270°F.  Figures  All  and  A12 

show  a  two  stage  cure  with  the  mold  set  at  300°F  (148. 9°C),  then  cooling  at 
room  temperature.  The  pad  Is  held  In  the  mold  long  enough  to  achieve  optimum 
cure  at  the  center.  Figure  A12  shows  that  this  results  In  severe  overcure 
near  the  outside  of  the  pad. 

So  far,  all  of  the  curing  cycles  have  been  for  the  case  of  steel  wire 
reinforcement  (X  direction,  Vf  =  0.013).  For  comparison,  the  cycle  depicted 
in  Figures  A7  and  A8  was  rerun  without  the  wire  reinforcement.  The  results 
were  that  the  cure  state  on  the  outside  did  not  charge  noticeably,  but  the 
cure  state  at  the  center  was  lower  by  roughly  20%.  See  Figures  A13  and  A14. 

In  a  last  example,  the  same  cycle  was  applied  to  a  pad  reinforced  with  the 
higher  conductivity  Copperply"  wire.  This  brings  the  cure  state  at  the 
center  substantially  closer  to  optimal.  See  Figures  AT 5  and  A16. 

Figures  A-17  through  A20  depict  the  following  case:  The  rubber  Is  Injection 
molded,  and  enters  the  mold  at  a  uniform  temperature  (100°C).  There  Is  no 
reinforcement.  The  case  of  30  minutes  In  the  mold  at  290°F  (143. 3°C), 
followed  by  60  minutes  of  post-heat  out  of  the  mold  at  100°C,  Is  shown  In 
Figures  A-17  and  A-18.  In  Figures  A-19  and  A-20,  the  mold  condition  Is  60 
minutes  at  133. 3°C,  10°C  less  than  the  previous  case.  Mold  time  was  doubled 
because  a  10°C  drop  In  mold  temperature  cuts  the  cure  rate  In  half  (see 
equation  34).  The  results  clearly  show  the  gain  In  cure  uniformity  resulting 
from  a  longer,  cooler  cure.  Figures  A- 21  and  A-22  show  a  variation  of  the 
second  case,  with  the  100°C  post-heat  dropped.  As  shown  before  for  the 
reinforced  rubber,  a  post-heat  at  a  temperature  lower  than  the  mold 
temperature  can  add  significantly  to  the  cure. 


,  .  .h;i+  rsn  he  draWn  from  these  models  of  curing  are  as  follows: 

The  conclusions  that  can  be  arV*  u  *n  severe  undercure  and/or  overcure.  2) 
1)  Simple  two  stage  cure  can  result  »  severe^  ^  un1form1ty  of  cure.  3) 

Preheating  prior  to  Riding  can  gr  y  i  P  preheating,  for  uniformity  of 
Post-heating  appears  to  b*  ^Iflcantly  less  time  In  the  mold.  4)  Curing 

cure.  It  can.  however.  ,f Icantly  Improve  the 

TtSe^re  5)  The  effect  of  metallic  reinforcement  on  curing  can 
be  Important. 


Listings  of  the  computer  codes  used 
output,  are  Included  In  Appendix  B. 
to  generate  the  figures. 


these  calculations,  plus  abbreviated 
The  output  Includes  tables  of  values  used 
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TABLE  A-  1 
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TEMPERATURE  (DEGREES  CENTIGRADE) 


^JFmIrTcal  simulatIonTtemperature  vs.  time  at  pad  center 
CYCLIC  LOADING,  HEATING  RATE  =  3.5E+05  ERG/(SEC*(CM**3)) 
VARIOUS  CONDS.  OF  REINFORCEMENT  AND  SURFACE  HEAT  TRANSFER 
POINTS  OFF  SCALE  ARE  STEADY  STATE  SOLUTIONS 


Figure  A-1  Hysteresis  Heating 


TEMPERATURE  (DEGREES  CENTIGRADE) 


NUMERICAL  SIMULATION  ..  TEMPERATURE  VS.  TIME  AT  PAD  CENTER 
CYCLIC  LOADING,  HEATING  RATE  REDUCED  BY  25%  (REINFCD  PADS) 
VARIOUS  CONDS.  OF  REINFORCEMENT  AND  SURFACE  HEAT  TRANSFER 
POINTS  OFF  SCALE  ARE  STEADY  STATE  SOLUTIONS 


Figure  A- 2  Hysteresis  Heating 
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TEMPERATURE  (DEGREES  CENTIGRADE) 


NUMERICAL  SIMULATION  ..  PAD  TEMPERATURE 
FUNCTION  OF  TIME  AND  LOCATION 
TWO  STAGE  CURE  CYCLE 


TIME  (MINUTES) 


Figure  A-3  Two  stage  Cure  Cycle,  Temperature  Vs  Time 


STATE  OF  CURE  (MINUTES) 


NUMERICAL  SIMULATION  ..  STATE  OF  CURE 
FUNCTION  OF  TIME  AND  LOCATION 
TWO  STAGE  CURE  CYCLE 


0  50  100  150  200  250  300 

TIME  (MINUTES) 


Figure  A-4  Two  Stage  Cure  Cycle,  Cure  State  vs  Time 


NUMERICAL  SIMULATION  ..  PAD  TEMPERATURE 
FUNCTION  OF  TIME  AND  LOCATION 
THREE  STAGE  CURE  CYCLE 
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STATE  OF  CURE  (MINUTES) 


NUMERICAL  SIMULATION  ..  STATE  OF  CURE 
FUNCTION  OF  TIME  AND  LOCATION 
THREE  STAGE  CURE  CYCLE 


o  50  100  150  200  250  300 

TIME  (MINUTES) 

Figure  A-6  Three  Stage  Cure  Cycle,  Cure  State  vs  Time 
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TEMPERATURE  (DEGREES  CENTIGRADE) 


NUMERICAL  SIMULATION  ..  PAD  TEMPERATURE 
FUNCTION  OF  TIME  AND  LOCATION 
THREE  STAGE  CURE  CYCLE 


Figure  A- 7  Three  Stage  Cure  Cycle,  Temperature  vs  Time 


CURE  (MINUTES) 


NUMERICAL  SIMULATION  ..  STATE  OF  CURE 
FUNCTION  OF  TIME  AND  LOCATION 
THREE  STAGE  CURE  CYCLE 


Figure  A-8  Three  Stage  Cure  Cycle,  Cure  State  vs  Time 


NUMERICAL  SIMULATION  ..  PAD  TEMPERATURE 
FUNCTION  OF  TIME  AND  LOCATION 
FOUR  STAGE  CURE  CYCLE 


CURE  (MINUTES) 


NUMERICAL  SIMULATION  STATE  OF  CURE 
FUNCTION  OF  TIME  AND  LOCATION 
FOUR  STAGE  CURE  CYCLE 


Figure  A- 10  Four  Stage  Cure  Cycle,  Cure  State  vs  Time 


TEMPERATURE  (DEGREES  CENTIGRADE) 


NUMERICAL  SIMULATION  ..  PAD  TEMPERATURE 
FUNCTION  OF  TIME  AND  LOCATION 
TWO  STAGE  CURE  CYCLE 


TIME  (MINUTES) 


Figure  A- 11  Two  Stage  Cure  Cycle,  Temperature  vs  Time 


CURE  (MINUTES) 


NUMERICAL  SIMULATION  ..  STATE  OF  CURE 
FUNCTION  OF  TIME  AND  LOCATION 
TWO  STAGE  CURE  CYCLE 
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CURE  (MINUTES) 


TEMPERATURE  (DEGREES  CENTIGRADE) 


TEMPERATURE  (DEGREES  CENTIGRADE) 


NUMERICAL  SIMULATION  ..  PAD  TEMPERATURE 
FUNCTION  OF  TIME  AND  LOCATION 
THREE  STAGE  CURE  CYCLE 
RUBBER  ONLY  (NO  REINFORCEMENT) 


LEGEND 
□  =  CENTER  (1,1) 

O  =  FAR  EDGE  (NI.l) 

A  =  NEAR  EDGE  (l.NJ) 

+  =  CORNER  (NI.NJ) 
x  =  B.C.  (OVEN  OR  MOLD) 


-ID- 


TIME  (MINUTES) 

Figure  A- 17  Injection  Mold,  Three  Stage  Cure,  Rubber  Only,  Temperature  Vs  Time 


NUMERICAL  SIMULATION  ..  PAD  TEMPERATURE 
FUNCTION  OF  TIME  AND  LOCATION 
THREE  STAGE  CURE  CYCLE 
RUBBER  ONLY  (NO  REINFORCEMENT) 


LEGEND 
□  =  CENTER  (1.1) 

O  =  FAR  EDGE  (NI,1) 

A  =  NEAR  EDGE  (l.NJ) 

+  =  CORNER  (NI.NJ) 
x  =  B.C.  (OVEN  OR  MOLD) 


Ed  120- 


TIME  (MINUTES) 

Figure  A-19  Injection  Mold,  Three  Stage  Cure,  Longer  Mold  Time,  Lower  Mold  Temp 


NUMERICAL  SIMULATION  ..  STATE  OF  CURE 
FUNCTION  OF  TIME  AND  LOCATION 
THREE  STAGE  CURE  CYCLE 
RUBBER  ONLY  (NO  REINFORCEMENT) 


LEGEND 

□  =  CENTER  (1,1) 

O  =  FAR  EDGE  (Nl.l) 

A  =  NEAR  EDGE  (1,NJ) 
+  =  CORNER  (NI.NJ) 
x  =  TARGET  VALUE 


ea 


|  TIME  (MINUTES) 

gure  A-20  Injection  Mold,  Three  Stage  Cure,  Longer  Mold  Time,  Lower  Mold  Temperature 


TEMPERATURE  (DEGREES  CENTIGRADE) 


NUMERICAL  SIMULATION  ..  PAD  TEMPERATURE 
FUNCTION  OF  TIME  AND  LOCATION 
TWO  STAGE  CURE  CYCLE 
RUBBER  ONLY  (NO  REINFORCEMENT) 


LEGEND 
□  =  CENTER  (1,1) 

O  =  FAR  EDGE  (NI,1) 

A  =  NEAR  EDGE  (l.NJ) 

+  =  CORNER  (NI.NJ) 
x  =  B.C.  (OVEN  OR  MOLD) 


75  100 

TIME  (MINUTES) 


Figure  A-21  Injection  Mold,  Two  Stage  Cure,  Temperature  Vs  Time 


STATE  OF  CURE  (MINUTES) 


NUMERICAL  SIMULATION  ..  STATE  OF  CURE 
FUNCTION  OF  TIME  AND  LOCATION 
TWO  STAGE  CURE  CYCLE 
RUBBER  ONLY  (NO  REINFORCEMENT) 


i  i  n i - 1  i  r 

0  25  50  75  100  125  150 

TIME  (MINUTES) 


Figure  A-22  Injection  Mold,  Two  Stage  Cure,  Cure  State  Vs  Time 
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APPENDIX  B 


COMPUTER  PROGRAMS 


o  n 


«***  TSO  FOREGROUND  HARDCOPY  **** 
QSNAME=ALGK610 . PAD20 . FORT 


REAL  KX,KY,K,KKR 

DIMENSION  CX(7),CY17)»CK(4),CH(2) >COR(  7)»T(7),T0(7) 

DIMENSION  X(  7  ) » Y(  7  ) 

C  2-0  CALC  OF  TEMP  IN  TRACK  PAD  . . . 

C  -A<X<+A,  -B<Y<+B 
C  ALL  UNITS  C.6.S. 

C  REF .  CARSLAU  AND  JAEGER,  CONDUCTION  OF  HEAT  IN  SOLIDS, 

C  SECTION  5.5,  P  170. 

C  GIVES  STEADY  STATE  TEMPERATURE  IN  A  2-0  BLOCK  WITH  INTERNAL  HEAT 
C  GENERATION  AND  CONVECTIVE  HEAT  TRANSFER 
C  KX  AND  KY  APE  THERMAL  CONDUCTIVITY  IN  X  AH)  Y  DIRECTIONS 
C  HX  AND  HY  ARE  COEFFS  OF  SURFACE  HEAT  TRANSFER  ON  X  AND  Y  SURFACES 
C  R  =  (1/HX)  =  SURFACE  THERMAL  RESISTANCE 
C  AO  =  RATE  OF  HEATING  PER  UNIT  VOLUME 
CX( I >=0. 

CY( 1  )  =  0. 

CX( 2 ) =0 . 5 
CY ( 2  )  =  0. 

CX( 3)=1.0 
CY(  3)=0. 

CX(4)=0. 

CY( 4  )  =  0 . 5 
CX( 5 )=0 . 

CY( 5  )  =  1 .0 
CX( 6 ) -0 . 5 
CY(6)=1.0 
CX( 7  )  =  1 . 0 
CY( 7  )  =  1 . 0 
PI  =  3. 1415927 

C  SIZE  OF  TEST  SPECIMEN  CROSS  SECTION  ( A*6 ) 

A=5 . 715 
8=2.54 
EPS=0 . 10 
KY  =  3 . 2E+04 
HY  =  2 . 35E+04 
AO=3.50E*05 

REDUCE  AO  TO  ACCOUNT  FOR  STIFFENING  DUE  TO  REINFORCEMENT 
A0=2 .625E+05 
C 

CK( 1 >=1.0 
CM  2  )  =  2 . 9375 
CK( 3 ) =5 .84375 
CK<4)=18. 

CH( 1 )=1.0 
CH( 2  )  =  2 . 0 
DOll  ICXY=1,7 
XC ICXY )=CX( ICXY )*A 
Y < ICXY  ) =CY( ICXY )*B 
11  CONTINUE 
DOl  IH= 1 , Z 
HX-CH(IH)#HY 
D02  IK= 1 , 3 
C  RESET  CORRECTION  TERMS 
DOiO  J=l,7 
COR( J)=0.00 


00000010 

00000020 

00000030 

00000040 

00000050 

00000040 

00000061 

00000062 

00000063 

00000064 

00000065 

00000066 

00000067 

0000006B 

00000070 

00000080 

00000090 

00000100 

OOOOOilO 

00000120 

00000130 

00000140 

00000150 

00000160 

00000170 

00000160 

00000190 

00000200 

00000210 

00000215 

00000220 

00000230 

00000240 

00000250 

00000260 

00000270 

00000271 

00000272 

00000273 

00000274 

00000280 

00000290 

00000300 

00000310 

00000320 

00000330 

00000340 

00000350 

00000360 

00000370 

00000380 

00000390 

00000400 

00000410 

00000420 

00000430 
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10  CONTINUE 

KX=CK(IK)*KY 
K=CK( IK )**0 .5 
R=1./HX 
RP=1./HY 
KKR=KX*RP/K 
CAL=i . /(R*KX ) 

WRITE* 6,120) 

120  FORMAT* IH1 ) 

WRITE*  6 1 1 0 1 ) 

101  FORMAT* //I OX, ' 2-0  THERMAL  MOOEL  OF  TRACK  PAD’//) 

WRITE*  6  » 102  )  A ,B » AO 

102  FORMAT*  I2X.  'A=*  ,E11  .4,5X>  1  B= 1  , Ell  .4 ,4X ,  1  A0=: 1  ,  E 1 1 . 4 ) 

WRITE* 6 ,10  3)  IH ,HX  »HY »R ,RP 

103  FORMAT* 2X, *IH=* ,12, 4X,,HX=I , El 1.4, 4X, *HY=*  ,E11.4, 

15X, ' Rs ' ,E11.4,4X, *  R  P=  *  » E 1 1 . 4 ) 

WRITE*  6,104)  IK,KX,KY,K 

104  FORMAT* 2X, 'IK=' ,I2,4X, 1 KX= 1 , El 1 .4 ,4X, ‘KY= * , E 1 1 .4 ,5X , *K= *  ,E1 1 .4 ) 
WRITE*  6,105)  KKR ,CAL 

105  FORMAT* 10X, * KKR=  * ,E11.4,3X, *CAL=  * ,Ell .4) 

C  N  IS  THE  IFCEX  OF  ALPHA 

003  N=l,7 
AN=N 

AL=l *2.*AN-1. )*PI/2 . -EPS )/A 
WRITE*  6,106) 

106  FORMAT* //10X, 'ITERATIVE  SOLUTION  FOR  AL*//) 

004  IA=1,8 

AAL~A*AL 

Z=AL*TAN(  AAD-CAL 

ZP=TAN* AAL )  +  AAlV  ( COS* AAL )**2 ) 

DEL =-Z/ZP 
AL=AL*0EL 

WRITE* 6 , 107 )  N,IA,Z,ZP, AL 

107  FORMAT*  3X,  ' N=  M2 , 2X,  •  IA=  M2 , 2X,  *Z=*  , E11.4,2X,  *ZP=  * » El  i  .4>2X> 

1 *  AL= ' ,E13,6> 

4  CONTINUE 
AAL=A*AL 

F 1  =  AL#AL*(  2 . *AAL+SIN*  2 . *AAL ) ) 

F2=KKR*AL*SIW*  K*B*AU*COSH*  K*B*AL 1 
F3=SIN*  AAL )/( F 1»F2  ) 

WRITE* 6 , 108 )  AALt F 1 , F2 »F3 

108  FORMAT*//10X, 'AAL=* ,Ell .4,2X, ^1=’ ,E11 .4,2X, *F2=’ ,Ell .4,2X, 
i’F3='  , El  1 .4 ) 

WRITE*  6,109) 

109  FORMAT* //I OX, 'SERIES  SOL  FOR  TEMP  AS  FN  OF  GEOMETRY'//) 

DOS  IC-l ,7 

TO* IC )=A0#< A*Rf *  A*A-X* IC )*X( IC ) )/*  Z . *KX ) ) 

COR* IC)=COR* ICH(4.*AO/KX)*F3*COS(X(IC)*AL)*COSH(K*Y( IC  )*AL) 
T*IC)=TO(IC)-COR(IC) 

WRITE* 6, 110)  N , IC ,X* IC ) , Y* IC ) ,T0( IC ) ,COR( IC  )  »T* IC  ) 

HO  FORMAT*  2X  ,  ' N=  '  ,  12 , 2X ,  *  IC-  1  , I2»2X>  *X= '  >F8.3,2X,'Y='  , F8 . 3, 2X ,  '  T0=  *  , 
1F8 . 3 , 2X ,  ' COR=  1  ,F8.3,2X,  ’T=SF8.3) 

5  CONTINUE 
3  CONTINUE 
2  CONTINUE 
1  CONTINUE 

STOP 

END 
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2-0  THERMAL  MODEL  OF  TRACK  PAO 


A*  0.5715E.01  B-  0.2540E+01  A0=  0.2625E+06 

IH=  l  HX=  0 . 2350E+05  HY=  0.2350E.05  R  =  0.4Z55E-04  RP=  0.4255E-04 

IK=  1  KX=  0 . 3200E.05  KY=  0.3200E+05  K=  0.1000E+01 

KKR=  0 .  1362E+01  CAL=  0.7344E+00 


ITERATIVE  SOLUTION  FOR  AL 


N= 

1 

IA= 

1 

Z  =  0 . 1831E+01 

ZP= 

0 . 1575E+03 

AL= 

0.245737E+00 

N= 

1 

IA= 

2 

Z=  0 . 7286E+00 

ZP= 

0.5714E+02 

AL= 

0 . 232984E+00 

N= 

1 

IA= 

3 

Z  =  0 . 2206E+00 

ZP= 

0 . 2780E+02 

AL= 

0 . 225049E+00 

N  = 

1 

IA= 

4 

Z=  0 . 3479E-0 1 

ZP= 

0 . 1973E+02 

AL= 

0 . 223285E+00 

M= 

1 

IA= 

5 

Z  =  0-1174E-02 

ZP  = 

0 . 1842E+02 

AL= 

0.223221E+0Q 

H= 

1 

IA= 

6 

Z=  0 . 1431E-05 

ZP= 

0  * 1837E+02 

AL= 

0 . 223221E+0Q 

N= 

1 

IA= 

7 

Z=  0 . 1Q13E-05 

ZP= 

0 , 1837E+02 

At  = 

0 . 223221E+0Q 

H= 

1 

IA  = 

8 

Z--0 . 1490E-05 

ZP= 

0 . 1837E+02 

AL= 

0 . 22322 1 E+OQ 

AAL=  0.1276E+01  Fl=  0.1549E+00  F2=  0.1347E+01  F3=  0.4567E+01 

SERIES  SOL  FOR  TEMP  AS  FN  OF  GEOMETRY 


N= 

1 

IC= 

1 

X= 

0.0 

Y= 

0.0 

10  = 

197.800 

COP  = 

150.523 

T= 

47.277 

N= 

1 

IC= 

2 

x= 

2.8S8 

Y- 

0.0 

TQ  = 

164.309 

COP  = 

120.926 

T= 

43.383 

N= 

1 

IC= 

3 

x= 

5.715 

Y  = 

0.0 

T0  = 

63.838 

COP  = 

43.775 

T= 

20.062 

N= 

1 

IC= 

4 

x= 

0.0 

Y= 

1.270 

T0  = 

197.800 

COR  = 

156.612 

T= 

41.188 

N= 

1 

ic= 

5 

x= 

0.0 

Y= 

2.540 

T0  = 

197.800 

COS  = 

175.372 

T= 

22.428 

N= 

1 

IC  = 

8 

x= 

2.858 

Y= 

2.540 

T0  = 

164.309 

COR  = 

140.889 

T  = 

23.420 

Hz 

1 

IC= 

7 

x= 

5.715 

Y  = 

2.540 

T0  = 

63.838 

co«= 

51.002 

T  = 

12.836 

ITERATIVE  SOLUTION  FOR  AL 


N= 

2 

IA= 

l 

Z=  0 . 7309E+01 

ZP= 

0.4727E+03 

AL= 

0.791605E+00 

N= 

2 

IA= 

2 

Z=  0 . 3418E+01 

ZP= 

0.1343E+03 

AL= 

0 . 766 147E+Q0 

N= 

2 

IA= 

3 

Z=  0 . 1475E+0 1 

ZP= 

Q.4366E+02 

Al= 

0.732372E+00 

N= 

2 

IA= 

4 

Z=  0 . 5246E+00 

ZP= 

0 . 1827E+02 

AL= 

0.703664E+00 

N= 

2 

IA= 

5 

Z=  0. 1 166E+00 

ZP= 

0.11 1 1E+02 

AL= 

0.693174E+00 

N= 

2 

IA  = 

6 

Z=  0 . 8358E-02 

ZP  = 

0 . 958 l E*0 1 

AL= 

0 . 692302E+00 

H= 

2 

IA= 

7 

Z=  0.4691E-04 

ZP= 

0 . 9470E+01 

AL= 

0.692297E+00 

N= 

2 

IA= 

8 

Z=  0.0 

ZP= 

0 ,  9469E  +  01 

AL= 

0.692297E+00 

AAL=  0.3956E.01  Fl=  0.4271E+01  F2=  0.5642E+01  F3=-0.3020E-01 

SERIES  SOL  FOR  TEMP  AS  FN  OF  GEOMETRY 

N=  2  IC=  1  X=  0.0  Y=  0.0  T0=  197.600  COR=  149.532  ' T=  46.268 


N-  2  TC=  Z  X*  2.656  Y=  0.0  TO*  164.309  COM  121.319  T« 

M=  2  IC=  3  X=  5.715  Y=  0.0  TO*  63.63d  COR*  44.455  T* 

N-  2  IC*  4  X*  0.0  Y*  1.270  TO*  197.600  COR*  155.213  T= 

N=  2  IC*  5  X*  0.0  Y*  2.540  TO*  197.600  COR*  172.411  T= 

N=  2  IC*  6  X*  2.656  Y*  2.540  TO*  164.309  COR*  142.062  T= 

N=  2  IC*  7  X*  5.715  Y*  2.540  TO*  63.636  COR*  53.033  T* 


ITERATIVE  SOLUTION  FOR  AL 


N*  3 

IA* 

1 

2= 

0 . 1279E+02 

ZP* 

0 . 7879E+03 

AL* 

0. 134055E+01 

N=  3 

IA* 

2 

2= 

0.61 34E+01 

ZP= 

0 . 21 39E+03 

AL* 

0 . 1 31 187E  +  01 

N=  3 

IA* 

3 

2= 

0 .27866*01 

ZP= 

0 . 6418E  +  02 

AL* 

0 . 1 26845E  +  01 

N=  3 

IA* 

4 

2= 

0. 1 101E+01 

ZP= 

0 . 2387E+02 

AL* 

0.122233E+01 

N*  3 

IA* 

5 

2= 

0 . 3003E+00 

ZP= 

0 . 1284E+02 

AL* 

0.11 9894E+0 1 

N=  3 

IA* 

6 

2= 

Q.3201E-01 

ZP= 

0 , 1029E  +  02 

AL* 

0. 119583E  +  01 

N*  3 

IA* 

7 

2= 

0.4063E-03 

2P* 

0 . 1003E  +  02 

AL* 

0 . 1 19S79E  +  0 1 

N*  3 

IA* 

8 

Z=- 

'0 . 5186E-05 

ZP= 

0. 1QQ3E  +  02 

AL* 

0 . 1 19579E+0 1 

AAL*  0.6834E+01  Fl  =  0.2082E+02  F2  =  0.2738E+02  F3=  0.9179E-03 

SERIES  SOL  FOR  TEMP  AS  FN  OF  GEOMETRY 


N= 

3 

IC* 

1 

X* 

0.0 

r= 

0.0 

TO* 

197.800 

COR* 

149.562 

T* 

46.236 

N* 

3 

IC* 

2 

X* 

2.656 

Y= 

0.0 

TO* 

164.309 

COR* 

121.290 

T= 

43.019 

N* 

3 

IC* 

3 

X* 

5.715 

Y= 

0.0 

TO* 

63.838 

COR* 

44.481 

T* 

19.357 

N* 

3 

IC* 

4 

X* 

0.0 

Y= 

1.270 

TO* 

197.800 

COR* 

*155.285 

T= 

42.515 

N= 

3 

IC* 

5 

X* 

0.0 

Y* 

2.540 

TO* 

197.800 

COR* 

172.726 

T* 

25.074 

N* 

3 

IC* 

6 

X* 

2.858 

Y  = 

2,540 

TO* 

164.309 

COR* 

141.760 

T= 

22.550 

N* 

3 

IC* 

7 

X* 

5.715 

Y  = 

2.540 

TO* 

63.838 

COR* 

53.301 

T* 

10.537 

ITERATIVE  SOLUTION  FOR  AL 


N* 

4 

IA* 

1 

2= 

0.1827E+02 

ZP= 

0. 1103E+04 

AL* 

0 . 188993E  +  0 1 

N* 

4 

IA* 

2 

2= 

0.8852E+01 

ZP* 

0.2936E+03 

AL* 

0 . 185979E  +  0 1 

N= 

4 

IA* 

3 

2= 

0 .4106E  +  Q1 

ZP= 

Q.8522E+02 

AL* 

0 . 181 162E+0 1 

N* 

4 

IA* 

4 

2= 

0 . 1668E+01 

ZP* 

0 . 3020E  +  02 

AL* 

0 . 175573E+0 1 

N* 

4 

IA* 

5 

2= 

0 . 4906E+00 

ZP* 

0.1562E+02 

AL* 

0 . 172431E+Q 1 

N  = 

4 

IA* 

6 

2= 

0 . 5568E-01 

ZP* 

0.1238E+02 

AL* 

0.171981E+01 

N* 

4 

IA* 

7 

Z* 

0 . 7463E-03 

ZP* 

0 . 1205E+02 

AL* 

0 . 171975E+Q 1 

N* 

4 

IA* 

8 

2= 

0 . 3576E-06 

ZP* 

0 . 1205E+02 

AL* 

0. 171975E*Q1 

AAL*  0.9828E+01  Fl=  0.6027E+02  F2=  0.1318E+03  F3=~0 . 4943E-04 

SERIES  SOL  FOR  TEMP  AS  FN  OF  GEOMETRY 


N= 

4 

IC* 

1 

X* 

0.0 

Y* 

0.0 

TO* 

197.800 

COR* 

149.560 

T = 

48.239 

N= 

4 

IC* 

2 

X* 

2.856 

Y* 

0.0 

TO* 

164.309 

COR* 

121.269 

T= 

43.020 

N= 

4 

IC* 

3 

X* 

5.715 

Y= 

0.0 

TO* 

63.636 

COR* 

44.482 

T  = 

19.356 

N= 

4 

IC* 

4 

X* 

0.0 

Y= 

1.270 

TO* 

197.800 

COR* 

155.277 

T* 

42.522 

N= 

4 

IC* 

5 

X* 

0.0 

Y* 

2.540 

TO* 

197.800 

COR* 

172.662 

T = 

25.138 

42.996 

19.363 

42.567 

25.366 

22.247 

10.605 


B-4 


N*  4  IC* 

6 

X* 

raw 

Y= 

2.540 

TO* 

144.309 

COP* 

141.747 

T* 

22.542 

N-  4  IC= 

7 

x= 

5.715 

Y  = 

2.540 

TO* 

63.636 

COP* 

53.340 

T= 

10.478 

ITERATIVE  SOLUTION  FOR  AL 

N= 

5 

IA* 

1 

2= 

0 . 2374E+02 

ZP* 

0 . 1418E+04 

AL= 

0.243945E+QI 

N= 

5 

IA* 

2 

2= 

0 . 1 157E+02 

ZP= 

0 . 3736E+03 

AL* 

0 . 24Q849E+Q1 

N= 

5 

IA= 

3 

2= 

0.5427E+01 

ZP* 

0. 1064E+03 

AL* 

0 . 235749E+01 

N= 

5 

IA* 

4 

2= 

0 . 2278E*01 

ZP* 

0 .  36  74E*02 

AL* 

0 . 229550E+0 l 

N= 

5 

I  A* 

5 

2= 

0 .6807E+00 

ZP* 

0. 1872E*02 

AL* 

0 . 225914E+0 1 

N= 

5 

IA  = 

4 

'  2= 

0 . 7654E-0 1 

ZP* 

0. 1493E  +  02 

AL* 

0.225402E+01 

N= 

5 

IA= 

7 

2= 

0 . 91 72E-03 

ZP* 

0.14S8E+02 

AL* 

0 . 225395E+0 1 

N* 

5 

IA= 

8 

2= 

0 . 2980E-06 

ZP* 

0 . 1457E+02 

AL* 

0 . 225395E+0 1 

AAL* 

Q.1288E+02  FI* 

0.1339E+03  F2= 

0.6235E+03  F3=  0.3711E-05 

SERIES 

SOL  FOP  TEMP 

AS  FN  OF 

GEOMETRY 

N= 

5 

IC* 

1 

X* 

0.0 

Y* 

0.0 

TO* 

197.800 

COR* 

149.560 

T= 

48.239 

H* 

5 

IC* 

2 

X* 

2.858 

Y* 

0.0 

TO* 

164.309 

COR* 

121.289 

T* 

43.020 

N= 

5 

IC  = 

3 

X* 

5.715 

r= 

0.0 

TO* 

63.638 

COR* 

44.482 

T* 

19.356 

N* 

5 

IC  = 

4 

X* 

0.0 

Y  = 

1.270 

TO* 

197.600 

COR* 

155.278 

T* 

42.521 

H= 

5 

IC= 

5 

X* 

0.0 

Y* 

2.540 

TO* 

197.800 

COR* 

172.681 

T= 

25.119 

N* 

5 

IC  = 

6 

x= 

2.858 

Y* 

2.540 

TO* 

164.309 

COR* 

141.765 

T* 

22.544 

N= 

5 

IC  = 

7 

X* 

5.715 

Y* 

2.540 

TO* 

63.838 

COR* 

53.378 

T* 

10.440 

ITERATIVE  SOLUTION  FOR  AL 

N* 

6 

IA= 

1 

Z=  0 • 2921E+02 

ZP= 

0 . 1732E  +  04 

AL= 

0 . 298904E+0 1 

N* 

6 

IA* 

2 

Z*  0 . 1429E+Q2 

ZP= 

0 .45  36E  +  Q3 

AL* 

0.295754E+01 

N= 

6 

IA* 

3 

Z*  0.6748E+01 

ZP  = 

0 . 1276E  +  03 

AL* 

0.290466E+01 

K= 

6 

IA* 

4 

Z=  0 . 2648E+0 1 

ZP* 

0 .  4337E  +  02 

AL* 

0 . 283854E+0 1 

N= 

6 

IA* 

5 

Z*  0 .8693E+00 

ZP* 

0.2197E+02 

AL* 

0 . 279896E+0 1 

N= 

4 

IA* 

6 

Z=  0.9504E-01 

ZP* 

0 . 1770E  +  02 

AL* 

0 . 279359E  *0 1 

N= 

6 

IA* 

7 

Z=  0 . 9795E-03 

ZP* 

0 . 1733E+02 

AL* 

0.279353E+01 

N= 

4 

IA* 

8 

Z=-0.1258E'04 

ZP* 

0 . 1733E+02 

AL* 

0.279353E+01 

AAL* 

0 . 1597E+02  Fi* 

0 . 2530E+03  F2* 

0 . 2898E+04  F3=-0 . 3447E-06 

SERIES  ! 

SOL  FOR  TEMP 

AS  FN  OF 

GEOMETRY 

N= 

6 

IC  = 

1 

X  = 

0.0 

Y  = 

0.0 

TO* 

197.800 

COR* 

149.540 

T* 

48.239 

N= 

4 

IC  = 

2 

X* 

2.858 

Y  = 

0.0 

TO* 

164.309 

COR* 

121.289 

T* 

43.020 

N= 

6 

IC* 

3 

X* 

5.715 

Y= 

0.0 

TO* 

63.838 

COR* 

44.482 

T* 

19.356 

N= 

4 

IC= 

4 

X* 

0.0 

Y= 

1.270 

TO* 

197.800 

COR* 

155.276 

T* 

42.522 

N= 

6 

IC* 

5 

X* 

0.0 

Y* 

2.540 

TO* 

197.800 

COR* 

172,674 

T* 

25.126 

N* 

4 

IC* 

6 

X* 

2.858 

Y* 

2.540 

TO* 

164.309 

COR* 

141.766 

T* 

22.543 

N* 

6 

IC* 

7 

X* 

5.715 

Y* 

2.540 

TO* 

63.638 

COR* 

53.384 

T* 

10.454 

B-5 


ITERATIVE  SOLUTION  FOR  AL 


N=  7  IA=  1 
N=  7  IA=  2 
N=  7  IA*  3 
N=  7  IA=  4 
M=  7  IA=  5 
N=  7  IA=  8 
M=  7  IA=  7 
N=  7  IA=  6 


IA=  1  Z=  0.3489E*02  ZP=  0.2048E*04 
IA=  2  Z=  0 . 1701E+02  ZP=  0.5338E+03 
IAs  3  Z=  0.8072E*01  ZP=  0.1489E+03 
IA=  4  Z=  0.3460E+01  ZP*  0.5006E+02 
;*  c  t  o  1058E+01  ZP=  0.2529E+02 
IA:  6  Z--  0  11212*00  ZP=  0 . 2058E+02 
A'  7  Z=  0.10172-02  Zps  0. 20212*02 
IA=  6  Z=-0. 82852-05  ZP=  0.20212*02 


AL=  0 . 353867E*01 
AU=  0 . 3506802*01 
AL=  0. 3452612*01 
AL=  0.3383502*01 
AL=  0 . 334167E+01 
AL=  0. 3338222*01 
AL=  0. 333617E*01 
Al=  0.333617E+01 


AAL=  0.1907E 


,02  Fl  =  0.4291E*03  F2=  0.1327E*05  F3=  0.3776E-07 


SERIES  SOL  FOR  TEMP  AS  FN  OF  6E0METRY 


N=  7  Its  l  X= 
N=  7  rc=  2  X= 
N=  7  IC=  3  X= 
H-  7  XC=  4  X= 
N=  7  IC=  5  X= 
H=  7  IC=  6  X= 
7  1C*  7  X= 


0.0 

2.858 

5.715 

0.0 

0.0 

2.858 

5.715 


0.0 

0.0 

0.0 

1.270 

2.540 

2.540 

2.540 


TQ=  197.800  COR- 
TO-  L64.309  COR  = 
T0=  63.838  COR- 

T0=  197.800  COR  = 
TO-  197.800  COR- 


COR-  149.560  T= 
COR-  121*289  T= 
COR=  44.482  T= 
COR-  155.278  T= 
COR-  172.877 


T0=  164.309  COR-  141.763  T= 
TO-  63.838  COR-  53.387  T= 


48.239 

43.020 

19.356 

42.521 

25.U3 

22.546 

10.451 


B-6 


I 


»***  T30  FOGEGROWC  HARDCOPY  **»* 
DSHAHE=AL6K610 . PA0CH20 . FORT 


C  2-0  MODEL  OF  A  RECTANGULAR  RL»BER  BLOCK,  USING  SYMMETRY 
C  MODIFIED  FOR  GRAPHICS  .... 

C  CAN  MOOEL  CURING  AMD  HYSTERESIS  HEATING 
C  CAN  USE  TEMPERATURE  AND  SURFACE  HEAT  TRANSFER  8.C. 

C  USES  FINITE  DIFFERENCE  METHOO  —  REF  WORK  OF  2/14/85 
C  UNITS  ARE  GRAM-CENTIMETER-SECOND 
C 

C  REFERENCES: 

C  HILLS »  "HEAT  TRANSFER  Ah©  VULCANIZATION  OF  RUBBER'1 
C  CARSLAW  AND  JAEGER,  "CONDUCTION  OF  HEAT  IN  SOLIDS" 

C 

C  VARIABLE  NAMES: 

C  DEGC  =  DEGREES  CENTIGRADE 

C  K  =  THERMAL  CONDUCTIVITY ,  GRAM*CM*( $**( -3 ) )*( DEGC**( -1  )  ) 

C  RHO  =  DENSITY,  GRAM#(  CM**( -3  )) 

C  C  =  HEAT  CAPACITY,  ( CM**2 )*( S»*( -2  )  )*< DEGC**( - 1  )  ) ,  =  ERG/( GPAM*DEGC ) 
Cl:  INDEX  IN  THE  X  DIRECTION 
CJ=  INDEX  IN  THE  Y  DIRECTION 

C  01  =  THERMAL  DIFFUSIVITY  IN  X( I )  DIRECTION,  ( CM**2 )»( S**< -I ) ) 

C  DJ  =  THERMAL  DIFFUSIVITY  IN  Y(J)  DIRECTION 
C  D  =  K/( RHO*C ) ,  GENERAL  DEFINITION  OF  THERMAL  DIFFUSIVITY 
C  HX  =  HEAT  TRANSFER  COEFFICIENT  IN  X  DIRECTION, 

C  HY  =  H  IN  Y  DIRECTION,  GRAM*( S**( -3 )  )#( OEGC**( - 1 ) ) 

C  W  =  HEATING  RATE,  GRAMM CM*«< -1 ) )*(  S( #*-3 ) ) ,  =  ERG/(  S*CM**3) 

C  TC  =  TEMPERATURE  COEFFICIENT  OF  VULCANIZATION 
C  VREF  =  REFERENCE  TEMPERATURE  FOR  CURING  IhCEX,  DEGC 
C  VO  =  INITIAL  TEMPERATURE  OF  THE  BLOCK 
C  NI  =  NUMBER  OF  NODES  IN  X  DIRECTION 
C  NJ  =  NUMBER  OF  NOOES  IN  Y  DIRECTION 
C  IBC  =  BOUNDARY  CONDITION  INDICATOR  ON  SURFACE  I  =  NI 

C  JBC  =  BOUNDARY  CONDITION  Ih©ICATOR  ON  SURFACE  J  =  NJ 

C  T  =  TIME,  SECONDS 
C  TM  =  TIME  IN  A  STEP,  SEC 

C  1  INDICATES  SURFACE  TEMPERATURE  BC  SPECIFIED 
C  2  INDICATES  CONVECTIVE  HEAT  TRANSFER  BC 
C  VS  =  SURFACE  TEMPERATURE  FOR  BC  TYPE  1 

C  VS  -  FREE  STREAM  TEMPERATURE  FOR  BC  TYPE  2 

C  VSI  =  SURFACE  TEMPERATURE  AT  BEGINNING  OF  STEP  ( INITIAL  TEMP. > 

C  VSF  =  SURFACE  TEMPERATURE  AT  END  OF  STEP  (FINAL  TEMP.) 

C  VST  =  SURFACE  TEMPERATURE  AT  TIME  T 
C  VSX  =  SURFACE  TEMPERATURE  ON  A  "X=CONST"  SURFACE 

C  VSY  =  SURFACE  TEMPERATURE  ON  A  "Y=CONST"  SURFACE 

C  HERE,  FOR  NOW,  VSX  =  VSY 
— C  LX  =  HALF  LENGTH  IN  X  DIRECTION,  CM 

C  LY  =  HALF  LENGTH  IN  Y  DIRECTION,  CM 

C  DX,  OY  =  CELL  SIZE 
C  DT  =  TIME  STEP 

C  A  =  TIME  STEP  CALCULATION  FACTOR 
C  V(I,J)  =  TEMPERATURE  AT  NOOE  I,J 
C  SC(I,J>  =  CURING  INDEX  AT  NOOE  I,J 
C  02VDX2C I , J )  =  SECOND  DERIVATIVE  OF  V  H.R.T.  X 
_C  SAME  IDEA  FOR  Y 
-C  NPPL  =  NUMBER  OF  POINTS  PLOTTED 
C  SCTV  =  TARGET  VALUE  OF  SC  (CURE  NUMBER )  (UNITS  ARE  TIME  (MINUTES)) 
REAL  KX,KY » LX, LY 

COMMON/A/  KX*KY ,RHO,C »DI ,DJ »HX »HY »W»DT,DX,DY >OXS »DYS , TC 


00000100 

ooooouo 

00000200 

00000300 

00000400 

00000500 

00000600 

00000610 

00000620 

00000630 

00001000 

00001100 

00001200 

00001300 

00001400 

00001500 

00001600 

00001700 

00001800 

00001400 

00002000 

00002100 

00002200 

00002300 

00002400 

00002500 

00002610 

00002700 

00002800 

00002900 

.00003000 

00003010 

00003020 

00003100 

00003200 

00003300 

00003400 

00003500 

00003600 

00003612 

00003620 

00003630 

00003640 

00003700 

00003800 

00003900 

00004000 

00004010 

00004100 

00004200 

00004300 

00004400 

00004410 

00004500 

00004600 


I 


B-7 


nnooonnooo  nooo 


9  VpFF.T 

COMMON/B/  NI.NJ.NW.NPWI.HPPl 

CONHON/C/  V(  11,11  ),S^!  }il1l1^JT,non),vlHJUOO).VNINJ(  100) 

DIMENSION  XDATAI  1000)  .YDATAl  1000),ZDATAl  1000 

SCTV=27 .5 

NPPL=0 

KY=3.2E+04 


KX=6.4E+04 


SPECIAL  CASE  ..  NO  HIRES  ... 

KX-KY 

FROM  CWW20  CALCULATION.  STEEL  WIRE  ... 

K.X=1 .2E*05 

CORRECTION  ..  SIMPLE  VOLLUE  AVERAGE.  STEEL  HIRE  ... 
KX=9.4E*04 

SPECIAL  CASE  ..  30X  CONDUCTIVITY  COPPERPLY  ... 
KX=1.87E+05 

RHO=l.lS 

C=1.8£*07 

OI=KX/lRHCH«C) 

OJ=KY/(RHO*C) 

C  MOVING  AIR  ..  CHANGE  HX  Ah®  HY 
C  HX=2.35E+04 

c  HY=2 . 35E+04 

C  STILL  AIR  .  •  CHANGE  HX  Ah®  HY 
HX=1.70E*04 
HY=l . 70E+04 

C  FOR  FINNED  SURFACE  ... 

C  HX-4.70E+04 

|  FOR  STIFFENING  EFFECT  OF  NIRES  ... 

c  WS2.62SE+05 

C 

C  FOR  CURING  CYCLE  ... 

M=0.00 

C 

”  TC=2.0 

VREF=143.3 

V0=20. 

NI=6 

NJ=6 

T=0. 

LX=5.715 
LY=2.54 
*1  AI=NI-l 

AJ=NJ-1 
OX=LX/AI 
DY=LY/AJ 


00004610 

00004620 

00004630 

00004640 

00004650 

00004845 


00004846 
00004848 
0000484? 
00004948 
00004949 
00004950 
00004951 
00004955 
00004956 
00004957 
00004958 
00004959 
00004960 
00004970 
00004980 
00004990 
00004991 
00005000 
00005100 
00005200 
00005300 
00005402 
00005403 
00005404 
00005405 
00005406 
00005407 
00005408 
00005409 
00005410 
00005411 
00005412 
00005420 
00005430 
00005440 
00005450 
00005460 
00005470 
00005612 
00005617 
00005700 
.4  00005800 
00005810 
00006000 
00006100 
00006600 
00007000 
00007100 
00007200 
00007300 
00007400 
00007500 


I 
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* 


I 


0XS=DX**2 
OYS=OY »*2 
001  1=1 »NI 
002  J= 1 »NJ 
V(I,J)=V0 
SC( I,J)=0.0 
2  CONTINUE 
1  CONTINUE 
A=0 . 2 

OT=A/< ( OI/DXS )  «■( DJ/DYS  )  ) 

C  ANOTHER  WAY  ... 

OTX=A*OXS/DI 
DTY=A*OYS/DJ 
0T=AMIN1 ( DTX»DTY  ) 

WRITE( 6 >11)  DTX,DTY,DT 

11  FORMAT! *  1 V//20X* 'CALCULATED  DTX,DTY>DT  = ' > 3( El  1 . 4 , 2X )// ) 
C 

DT=25. 


C 

WRITE! 6 1 12  )  OT 

12  FORMAT( 20X> 'MOOIFIED  DT  =’»E11.4//) 

WRITE! 6  » 14  )  OX »0Y 

14  FORMAT( 20X> ' CELL  SIZE:  OX  = ' , El 1 .4 , 3X, ' OY  =',E11.4//) 
WRITE! 6  » 13 )  DI »0J 

13  FORMAT(  20Xt  ’  OIFFUSIVITIES :  DI  =  ' , El  1 .4 , 3X , ' 0 J  =’,E11.4//) 
NPR  =  0 

NPRM=12 

MT=3 

TfIC  1  )=7200 . 

VSI ( 1 )=100. 

VSF( 1 )=100. 

IBC( 1 )=2 
JBC( 1  )  =  2 
TM( 2 )=3900 . 

VSI( 2 )=132.22 
VSF( 2 )=132.22 
IBC!2)=i 
JBC! 2  )  =  1 
TM(  3 )=4200 . 

VSI! 3 )=20. 

VSF ( 3  )  =  20 . 

IBC( 3  )=2 
JBC( 3 )  =  2 
TM( 4 ) =1800 . 

VSI( 4  )  =  20 . 0 
VSF ( 4  )  =  20 « 0 
IBC!  4  )=2 
JBC(  4  )=2 
TT1!5)=3600. 

VSK  5  I  — 143 . 3 
VSF ( 5 ) =20 . 0 
IBC(  5  )=2 
JBC( 5  )=2 
MT=2 

TM< 1 )=4500. 

VSK  1  )=148.69 
VSF( 1 )=146,89 
IBC( 1  )  =  1 
JBC( 1 )=1 
TM(2)=7200. 


00007600 

00007700 

00008000 

00008100 

00008200 

00008300 

00008400 

00008500 

00008501 

00008502 

-00008503 

00008504 

00008505 

00008506 

00008507 

00008508 

00008509 

00008510 

00008511 

00008512 

00008513 


00008514 

00008515 

00008516 

00008517 

00008532 

00008534/- 

00008535 

00008536 

00008537 

00008538 

00008539 

00008540 

00008541 

00008542 

00008543 

00008544 

00008545 

00008546 

00008547 

00008548 

00003549 

00008550 

00008551 

00008552 

00008553 

00008554 

00008555 

00008556 

00008557 

00008558 
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VSI<2)=20. 
VSF(2)=20. 
IBC12)=2 
JBC( 2  )=2 


MT=  1 

TMC 1)=7200. 

VSItl)=20. 

VSFt i)=20. 

IBC(1)=2 
JBC(1)=2 
WRITEt 6,101 ) 

101  FQRMAT(  '1V/10X,  ******  2-0  BLOCK  CYCLE  *****’//) 

WRITEt  6  » 105 )  LX,  LY 

105  FORttATt  10X,  *  LX,  LY  = ' ,2t El 1 .4, 2X >/ ) 

WRITEt  6  » 102 )  OXiOY »DT 

102  FORttATt  10Xi  '  DX»  DY,  DT  s* ,3( El  1 .4,2X  )// ) 

WRITEt  6  » 106 )  KXtKY 

106  FORttATt 10X, *KX,  KY  = 1 , 2( Ell ,4 ,2X )/ ) 

WRITEt  6, 104 )  01, OJ 

104  FORttATt 10X, '01,  OJ  = ' , 2( El 1 .4 , 2X )/ ) 

WRITEt  6 » 103  )  HX.HY 

103  FORttATt 1 OX, 'HX,  HY  = * ,4( Ell .4»2X )//  ) 

WRITEt  6, 111) 

111  FORttATt //I OX, ' TIME -TEMPER ATURE-B.C. -HEATING  HISTORY1//) 

0022  NT=1,MT 

WRITEt 6, 112)  NT ,TM(NT) ,VSItNT ) ,VSF(  NT ) , IBCt NT ), JBCt  NT  )  ,M 

112  FORttATt  5X,  'NT=*  »I2»2X»  * TTi( NT )  =  '  ,F7.1,2X,  *  VSIt  KT )=  •  , F7. 1 ,2X 
1 , 1  VSF(MT)= 1  ,F7. 1,2X, 1  IBCtMT )s ’  ,11  »2X, 1  JBCtNT )  = '  ,11 , 2X,  'W=* 

2, Ell. 4/) 

22  cotmwE 
003  NT=l ,MT 

CALL  STEPt  TMtKT ) , VSIt NT ) ,VSFt  NT ), IBCt  NT  ) , JBCt  NT  )  ) 

3  CONTINUE 
WRITE(6 ,71 )  NPPL 

71  FORMATt  *  1 V/10X, '***  PLOTTING  OUTPUT,  NPPL  =*,I3//) 

WRITEt  6 ,72 ) 

72  FORttATt  16X,  ‘HP*  ,3X,  *T( MINUTES)*  ,6X,  ‘  VtB.C.  ) '  ,8X,  '  Vll/'  ,8X*  *VNI1/' 
2,8X, ' V1NJ/' ,7X, 'VNINJ/1 ) 

WRITEt  6,73) 

73  FORttATt  5 IX, *  SCI 1 1 ,8X, 'SCNIl* ,8X, *SC1NJ’ , 7X, ’ SCNINJ * //  ) 

0074  NP=1,NPPL 

WRITEt  6, 75)  NP,TP(NP),VBC(NP),V11(NP) , VNI1 ( NP ) , VINJt NP ) ,  VNINJ  tNP) 

75  FORttATt 15X, 13, 6(2X, Ell. 4)) 

WRITEt 6, 76)  SCI 1(NP), SCNIl t NP ), SCINJt NP ) .SCNINJ(NP) 

76  FORttATt 44X, 4 ( 2X,Ei 1 .4 )/ ) 

74  CONTINUE 
NEX=G 

00201  NL=1,5 
00202  NP=1,NPPL 
NEX=NEX+1 
XDATAt NEX )=TP( NP ) 

IFtNL.EQ. 1 )  YOATAt  NEX )=V11 ( NP ) 

IF( NL.EQ.2 )  YOATAt  NEX )=VNI1 ( NP ) 

IFtNL.EQ. 3)  YOATAt  NEX )=V1NJ( NP ) 

IFtNL.EQ. 4)  YOATAt  NEX )=VNINJ( NP ) 

IFtNL.EQ. 5)  YOATAt NEX )=VBCtNP) 

IFtNL.EQ. 1)  ZDATAt  NEX  )=SC1 it  NP ) 

IFtNL.EQ. 2)  ZDATAt NEX )=SCNI 1 ( NP ) 

IFtNL.EQ. 3)  ZDATAt NEX )=SClNJtNP) 

IFtNL.EQ. 4)  ZDATAt  NEX )=SCNINJ( NP ) 


00008562 

00008563 


00008564 

00008570 


00008571 
00008572 
00008580 
00008590 
00008591 
00008592 
00008593 ' 
00008594 
00008595 
00008596 
00008597 
00003598 
00008599 
00008600 
00008610 
00008620 
00008630 
00008631 
00003640 
00008650 
00008660 
00008670 
00008680 
00008690 
00008691 
00008692 


I 


ft 


IF(NL . EQ.5 )  ZD AT A! HEX )=SCTY 

202  CONTINUE 
201  CONTINUE 

WRITE! 6  »203 ) 

203  FORMAT! *  1V/1QX,  ******  EZGRAF  INPUT  **#**•//) 

00204  NE  =  1 »NEX 

WRITE ( 6 , 205 )  NE , XD ATA! NE  ) , Y0ATA( NE ) , ZD ATA( NE  ) 

205  FORMAT! 10X> 'NE=' »I4,3X» 'XDATA='  »EU.4»3X» 'YDATA=' »E11.4»3X 
2  * ZDATA= ' i El 1 . 4 ) 

204  CONTINUE 
WRITE! 6*210) 

210  FORMAT! * 1 *//10X » ' FORMATTED  PARTIAL  EZGRAF  INPUT’//) 

WRITE! 6  » 206  ) 

206  FORMAT! ’  XDATA '  ) 

WRITE ( 6 » 207  )  (  XDATA! N ) ,N= 1 ,NEX ) 

207  FORMAT! 8X»E1 1 .4, 2X> Ell . 4,2X,E1 1 .4,2X»E1 1 . 4»2X»E11.4  ) 

WRITE ( 6  » 208  ) 

208  FORMAT! 1  YDATA* ) 

WRITE (6*207)  ( YDATA! N ) ,N= 1 ,NEX ) 

WRITE ( 6  *  209 ) 

209  FORMAT! •  ZD AT A*  ) 

WRITE! 6 * 207 )  ( ZDATA! N ) ,N=1 ,NEX ) 

STOP 

END 

SUBROUTINE  STEP! TM , VSI , VSF , IBC  » JBC ) 

REAL  KX,KY 

COMMON/ A/  KX»KY >RHO»C  >01 *DJ  >HX,HY ,W»OT»DX >DY  >DXS,DYS»TC 
2  » VREF  i  T 

COMMOtl/B/  NI  >NJ  ,NPR  »NPRM  ,NPPL 
COMMON/C/  V(lUU),SC(liai) 

COMMON/D/  TP!  100  )  ,VU(  100 )  ,VNI1(  100 )  ►  V1NJ!  100  )  ,VNINJ(  100) 
COMMON/E/  SCll! 100)>SCNIi< 100),SC1NJ( 100),SCNINJ! 100)»VBC( 100) 
DIMENSION  02VDX2! ll»ll)  ,D2VDY2! 11,11) »DSCOT( 11,11) 

NTSUB=TM/DT 
TSUB=0.Q 
DELVS=VSF-VSI 
DO!  NT=  1 » NT  SUB 

VS=VSI +0ELVS*! TSUB+0 . 5*0T )/TM 

vsx=vs 

VSY=VS 
0011  1=1, NI 
0012  J= 1 ,NJ 

IF(I.EQ.l)  D2VDX2! I»J)=2.*!V(2»J)-V(  1 ,  J  ) )/DXS 
IF!(I.NE.i).AND.(I.NE.Nin  02VDX2! I, J )  =  ( V!  ( I ♦  1 ) ,  J  )-2,*V!  I 
2»J)+V((I-1),J)  )/DXS 

IF! ( I.EQ.NI ) .AM). ! IBC.EQ.2 ) )  VOUT=V! (NI-1 ) • J )-2 .*DX*HX*( V 
2! NI , J ) -VS )/KX 

IF! (I.EQ.NI). AND.! IBC.EQ.2))  02V0X2(NI * J )=! VOUT-2 . *V( NI » J 
2)+V((NI-l),J) )/DXS 

IF! (I.EQ.NI).  AND.  (IBC.EQ.m  D2VDX2! NI , J )=9999. 

IF! J.EQ.l)  D2VDY2! I»1  )=2.*( V(I,2)-V!  I>1 ) )/DYS 
IF!  C  J.NE.U.AKO.!  J.NE.NJ))  02VDY2!  I,J)  =  !V(I,(JM)  )-2.*V(I 
2»JNV!I,(  J-l)))/OYS 

IF!  ( J.EQ.NJ)  .AND.  (  JBC.EQ.2 ) )  VOUT=V(  I,  ( NJ-1 )  )-2 . *DY*HY*(  V 
2! I »NJ )-VS )/KY 

IF!  (  J.EQ.NJ ). AND. (JBC.EQ.2))  02VDT2! I»NJ )=( VOUT-2 . *V( I ,NJ 
2  )*V( I » ( NJ-1 ) ) )/DYS 

IF! (J.EQ.NJ ). AND. (JBC. EQ.l ) )  02VDY2! I »NJ  >=9999. 

12  CONTINUE 
11  CONTINUE 


00008700 

00008800 

00008900 

00009000 

00009010 

00009020 

00009310 

00009320 

00009330 

00009340 

00009400 

00009500 

00009600 

00009700 

00009800 

00009900 

00010000 

00010100 

00010191 

00010192 

00010193 

00010194 

00010195 

00010196 

00010197 

00010198 

00010199 

00010200 

00010201 

00010202 

00010203 

00010204 

00010205 

00010206 

00010207 

00010208 

00010209 

00010210 


I 


I 
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0021  1=1 *NI 

0V=fol*O2VDX2t I ,  J  )+DJ»D2VDY2t I  *  J )  + 1 W/l RHO*C ) ) )*0T 
va,ji=va.J)*ov 

22  CONTINUE 
21  CONTINUE 

IF(IBC.NE.l)  GO  TO  24 
0023  J=1.NJ 
VINI > J )=VSX 

23  CONTINUE 

24  IF(JBC.NE.l)  GO  TO  26 
0025  1=1 >NI 

V( I.NJ >=VSY 

25  CONTINUE 

26  CONTINUE 
005-1  1  =  1. NI 
0052  J=1  >NJ 

DSCOT(I.J)=«TC»hM0.1»<V!I.J)-VREF)>> 

NOTE  ..  OT/60 .  SO  SC(I.J)  HAS  UNITS  OF  MINUTES  ... 

SCI  I. J )=SCl I > J )+OSCDTt I, J )*DT/60 . 

52  CONTINUE 
51  CONTINUE 

TSUB=TSUB*DT 

VST=VSI+OELVS»TSUe/TM 

T=T+DT 

NPR=NPR+1 

IF1NPR.NE.NPRH)  GO  TO  41 

WRITE! 6 >31 1  T  ,  ro  .... 

31  FORMAT! ' 1 ' //50X, 'TIME! SEC )=  >F6.2//) 

WRITE!  6  >36  )  VST 

36  FORMAT! 50X, ’VS  AT  TIME  T  IDES.  Cl-  .F7.2//1 

32  FORMATISOxI'BLOCK  TEMPERATURE  (DEGREES  C)'//) 

33  FORMAT?  12x! '  1  =  1  ‘  >7X.  ’  1=2’  >7X, ' >7X>^I=4'  ’Z*) 

Z ' 1=6’ >7X, ’ 1=7' >7X> ' 1=8' >7X» ' 1=4  >6X,  1-10  .6X,  1-11  /> 

0027  J=1 >NJ 

WRITE(6,34)  J. !  V!  I .  J )  >1=1  >NI ) 

34  FORMAT! 3X,'J=’ ,12 ,2X, 11! F6.2 .2X) ) 

27  CONTINUE 

WRITE! 6 >6 1  1  ,  „ 

61  FORMAT! //50X,  ’RATE  OF  CURE'//) 

WRITE!  6.33) 

0062  J=1 ,NJ  .  .  M_, 

WRITE! 6 >63 )  J, ! OSCOT! I , J 1 .1  =  1 .NI 1 
63  FORMAT! 3X> ’ J= ‘ . 12 >3X > 111 F7.5.3X1/1 

62  CONTINUE 
WRITE16  >35 1 

35  FORMAT! ///50X, ’CURE  STATE  (MINUTES)  //) 

WRITE! 6. 33) 

0026  J=l.NJ 

WRITEI6.37)  J>!SC( I > J  )  >1=1 >NI) 

37  FORMAT! 3X, 1 J=’ ,I2»3X. 11! F8.3.2X) ) 

26  CONTINUE 
I  NPPL=NPPL*1 

TP! NPPL )=T/60 . 

VBC(NPPU=VS 
VU!NPPL)=V(  1,1) 


00010220 

00010221 

00010222 

00010223 

00010226 

00010227 

00010226 

00010229 

00010230 

00Q10231 

00010232 

00010233 

00010234 

00010235 

00010236 

00010237 

00010236 

00010239 


00010240 

00010241 

00010242 

00010243 

00010244 

00010245 

00010246 

00010250 

00010310 

00010320 

00010331 

00010332 

00010333 

00010340 

00010350 

00010360 

00010370 

00010360 

00010390 

00010391 

00010392 

00010393 

00010394 

00010395 

00010396 

00010397 

00010398 

00010400 

00010401 

00010402 

00010403 

00010404 

00010405 

00010406 

00010407 

00010408 

00010409 


K- 


I 
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VNI1<NPPL)=V(NI,1> 
V1NJ(NPPL)=V(1,NJ> 
VNINJ(NPPl>=V(NI,NJ) 
SC11(NPPU=SC(  1.1) 
SCNIltNPPL )=SC(  HI > 1 ) 

SC1NJ(NPPL)=5C( l.HJ) 

SCNINJl  HPPD=SC(  NI  *NJ  ) 
CONTIHUE 

IFINPR.EQ.NPRM)  NPR  =  0 
CONTINUE  . 

RETURN 

END 


00010410 

00010411 

00010412 

00010413 

00010414 

00010415 

000-10416 

00010417 

00010410 

00010420 

00010500 

00010600 


CALCULATED  DTX.DTY >0T  =  0.1690Et03 


0 . 3339E+02  0.3339E+02 


MOOIFIED  DT  =  0.2500E+02 

CELL  SIZE:  DX  =  0.1143E+01  DY  =  0.5080E+00 
DIFFUSIVITIES:  DI  *  0.1546E-02  DJ  =  0.1546E-02 


I' 


I 
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*****  2-0  BLOCK  CYCLE  ***** 


LX.  LY  =  0.5715E+01  0.2540E+01 

OX,  OY,  DT  =  0.I143E+01  0.5080E+00  0.2500E*02 

KX,  KY  =  0.3200E+05  0.3200E+0S 

01,  OJ  =  0.1546E-02  0.1546E-02 

HX,  HY  =  0.1700E+05  '  0.1700E+05 


TIME-TEMPERATURE-B.C. -HEATING  HISTORY 


NT= 

l 

Tf1(  NT  )  =  7200.0 

VSI(NT)  = 

100.0 

VSF ( NT  )  = 

100.0 

IBC(  KT  1=2 

JBCtNT  )=2 

14= 

0.0 

HI¬ 

2 

TM(NT)=  3900.0 

VSI(  NT  )  = 

132.2 

VSF(NT)= 

132.2 

IBCCNT)  =  1 

JBC(NT)  =  1 

Ms 

0.0 

NTS 

3 

TH(NT)=  4200.0 

VSI(NT)  = 

20.0 

VSF<NT)= 

20.0 

IBCt  NT  )=2 

JBC(KT)=2 

0.0 
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TIME(SEC)=  300.00 
VS  AT  TIME  T  IDES.  C)=  100.00 
BLOCK  TEMPERATURE  (DEGREES  C) 


1=1 

20.17 

20.42 

21.57 

24.82 

31.99 

44.48 

1=2 

20.  IT 
20.42 
21.57 
24.62 
31.99 
44.46 

1=3 

20.19 

20.44 

21.59 

24.64 

32.01 

44.49 

1=4 

20.51 

20.76 

21.91 

25.15 

32.29 

44.73 

1=5 

23,61 

23.85 

24.96 
28.10 

34.97 
46.92 

1=6 

41.68 

42.07 

42.94 

45.35 

50.61 

59.70 

RATE  OF 

CURE 

1=1 

0.00020 

0.00020 

0.00022 

0.00027 

0.00045 

0.00106 

1=2 

0.00020 

0.00020 

0.00022 

0.00027 

0.00045 

0.00106 

1=3 

0.00020 

0.00020 

0.00022 

0.00027 

0.00045 

0.00106 

1=4 

0.00020 

0.00020 

0.00022 

0.00026 

0.00046 

0.00108 

1=5 

0.00025 

0.00025 

0.00027 

0.00034 

0.00055 

0.00126 

1=6 

0.00068 

0.00090 

0.00095 

0.00113 

0.00162 

0.00304 

CURE  STATE  (MINUTES) 


TIME ( SEC  )  =  600.00 


1=1 

22.18 

23.04 

25.85 

31.17 

39.59 

51.31 


1=1 

0.00023 

0.00024 

0.00029 

0.00042 

0.00076 

0.00170 


VS  AT  TIME  T  (DES.  C)  =  100.00 


BLOCK  TEMPERATURE  (DEGREES  C) 


1=2 

1  =  3 

1=4 

1-5 

1=6 

1  =  7 

22.20 

23.07 

25.88 

31.20 

39.62 

51.33 

22.45 

23.31 

26.11 

31.42 

39.81 

51.48 

23.96 

24.81 

27.56 

32.76 

40.99 

52.44 

30.  64 
31.61 
34.12 
36.87 
46.36 
56.77 

51.67 

52.41 

54.16 

57.47 

62.69 

69.93 

RATE  OF 

CURE 

1=2 

1=3 

1=4 

1=5 

1=6 

1=7 

0.00023 

0.00024 

0.00029 

0.00042 

0.00076 

0.00170 

0.00023 

0.00024 

0.00030 

0.00043 

0.00077 

0.00172 

0.00026 

0.00027 

0.00033 

0.00047 

0.00083 

0.00184 

0.00041 

0.00043 

0.00052 

0.00072 

Q.QOUl 

0.00246 

0.60177 

0.00164 

0.00207 

0.00261 

0.00374 

0.00619 

CURE  STATE  (MINUTES) 


TIMEC  SEC  »=  7200.00 


1=1 

64.02 

66.27 

67.03 

66.26 

89.92 

91.94 


1=1 

0.01866 

0.01920 

0.02023 

0.02203 

0.02472 

0.02845 


VS  AT  TIME  T  (DES.  C>=  100.00 


1=2  1=3  I=« 

86.40  87.56  89.42 
86.65  87.78  69.61 
87.39  86.46  90.18 
86.58  89.55  91.11 
90.20  91.03  92.37 
92.17  92.83  93.90 


BLOCK  TEMPERATURE  (DEGREES  C) 


1  =  5 

1=6 

91.90 

94.86 

92.05 

94.96 

92.48 

95.23 

93.20 

95.69 

94.16 

96.30 

95.33 

97.04 

RATE  OF 

CURE 

1=2 

0.01938 

0.01971 

0.02074 

0.02254 

0.02520 

0.02889 


1=3 

0.02099 

0.02132 

0.02234 

0.02410 

0.02669 

0.03025 


1=4 

0.02388 

0.02420 

0.02518 

0.02686 

0.02930 

0.03256 


1  =  5 

0.02835 

0.02865 

0.02953 

0.03103 

0.03316 

0.03598 


1=6 

0.03482 

0.03505 

0.03573 

0.03687 

0.03846 

0.04050 


1  =  7 


CURE  STATE  (MINUTES) 


TXME( SEC  )  = 

11100.00 

VS  AT  TIME 

T  (OES.  C )=  132.22 

BLOCK  TEMPERATURE  (DEGREES  C) 

1=1 

127.81 

128.03 

128.65 

129.63 

130.86 

132.22 

1=2 

128.01 

128.22 

128.82 

129.75 

130.92 

132.22 

1  =  3 

128.61 

128.79 

129.30 

130.10 

131.10 
132.22 

1=4 

129.57 

129.70 

130.08 

130.66 

131.40 

132.22 

1=5 

130.81 

130.88 

131.08 

131.39 

131.79 

132.22 

1=6  1=7 

132.22 

132.22 

132.22 

132.22 

132.22 

132.22 

BATE  OF  CU?E 

1=1 

0.34180 

0.34695 

0.36234 

0.38768 

0.42214 

0.46394 

1=2 

0.34681 

0.35160 

0.36646 

0.39088 

0.42397 

0.46394 

1=3 

0.36126 

0.36571 

0.37894 

0.40050 

0.42943 

0.46394 

1=4 

0.38609 

0.38958 

0.39988 

0.41646 

0.43834 

0.46394 

1=5 

0.42088 

0.42290 

0.42879 

0.43813 

0.45018 

0.46394 

1=6  1-7 

0.46394 

0.46394 

0.46394 

0.46394 

0.46396 

0.46394 

CURE  STATE  (MIHJTES) 

1=1 

10.899 

11.326 

12.722 

15.498 

20.693 

31.623 

1=2 

11.194 

11.621 

13.011 

15.766 

20.899 

31.658 

1  =  3 

12.200 

12.622 

13.990 

16.672 

21.595 

31.773 

1=4 

14.358 

14.767 

16.077 

18.588 

23.055 

31.999 

1=5 

18.922 

19.280 

20.404 

22.476 

25.947 

32.406 

1=6  1-7 

32.226 

32.257 

32.351 

32.519 

32.777 

33.157 

TIME ( SEC  )=15300.00 


4 


VS  AT  TIME  T  (DEG.  C)  =  20.00 


BLOCK  TEMPERATURE  (DEGREES  C> 


n 

M 

J= 

1 

64.81 

J= 

2 

64.00 

J= 

3 

61.57 

J= 

4 

57.63 

J= 

5 

52.31 

J= 

6 

45.82 

1  =  2  i  =  : 

63.76  60 

62 . 96  59 

60.59  57 

56.74  54 

51.55  49 

45,21  43 


;  i=4 

52  54.98 
78  54.34 
59  52.45 
02  49.37 
21  45.22 
34  40.15 


1  =  5 

1=6 

47.14 

37.30 

46,64 

36.98 

45.17 

36.05 

42.79 

34.52 

39.57 

32.47 

35.63 

29.97 

1=7  1=8 


1=9  1=10 


RATE  OF  CURE 


1  =  1 

1=2 

1  =  3 

J  = 

1 

0.00434 

0.00403 

0.00322 

J  = 

2 

0.00410 

0.00381 

0.00306 

J  = 

3 

0.00347 

0.00324 

0.00263 

J  = 

4 

0.00264 

0.00248 

0.00205 

J  = 

5 

0.00182 

0.00173 

0.00147 

j  = 

6 

0.00116 

0.00111 

0.00098 

1=4 

1  =  5 

1=6 

0.00219 

0.00127 

0.00064 

0.00210 

0.00123 

0.00063 

0.00184 

0.00111 

0.00059 

0.00149 

0.00094 

0.00053 

0.00112 

0.00075 

0.00046 

0.00079 

0.00057 

0.00039 

1=7  1=8  1=9  1=10 


1=1 


1  =  1 


CURE  STATE  (MIt*JTES) 


1=1 

1  =  2 

1  =  3 

J=  1 

18.294 

18.540 

19.278 

J=  2 

18.274 

18.526 

19.294 

J=  3 

18.468 

18.730 

19.552 

J=  4 

19.614 

19.868 

20.691 

J=  5 

23.091 

23.290 

23.948 

J=  6 

32.512 

32.543 

32.640 

1=4 

1  =  5 

1=6 

20.541 

22.980 

33.298 

20.638 

23.196 

33.312 

21.070 

23.893 

33.352 

22.290 

25.245 

33.410 

25.266 

27.723 

33.462 

32.815 

33.081 

33.469 

1=8  1=9  1=10  1=1 
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***  PLOTTING  output,  NPPL  =  51 


HP 

T(  MINUTES) 

V(  B .  C .  ) 

Vll/ 

SC11 

VNI 1/ 
SCNI1 

V1NJ/ 

SC1NJ 

VNINJ/ 

SCNINJ 

1 

0 . 5000E+0 l 

0. 1000E*03 

0.201 7E  +  02 

0 . 9733E-03 

0.4188E402 

0.2717E-02 

0 . 4448E  +  02 
0.3570E-02 

0.5970E402 

0.8199E-02 

2 

0 . 1000E+02 

0 . 1000E4-03 

0.2218E402 

0.2020E-02 

0.S187E*02 

0 . 9470E - 02 

0.5131E*02 

0 . 1 063E-0 1 

0.6993E+02 
0 . 3204E-0 1 

3 

0 . 1500E  +  02 

0.1000E403 

0 . 2583E+02 

0 . 331 1 E-02 

0 . 5862E+02 
0.21 10E-01 

0 .5570E+Q2 
0.2078E-01 

0 . 7530E+02 
0 . 7067E-0 1 

4 

0.2000E*02 

0 . 1000E  +  03 

0.3012E402 

0 . 5024E-02 

0 .6386E  +  02 

0 . 3850E-01 

0 .5910E+02 

0 . 3397E-01 

0. 78868*02 
0 . 1224E*00 

5 

0 . 2500E+02 

0 . 1 000E  +  03 

0.3457E402 

0.7348E-02 

0 . 6814E+02 

0 .  6254E-0 1 

0.62Q5E*02 

0.5033E-01 

0.8152E*02 
0 . 1862E*00 

4 

0 . 3000E+02 

0 . 1000E  +  03 

0 . 3698E+02 

0. 1051E-01 

0 . 7170E+02 

0 . 9394E-01 

0.6474E+02 

0 . 7019E-01 

0.8365E*02 
0.261 2E+00 

7 

0 . 3500E+02 

0 . 1 000E  +  03 

0.4325E402 

0.1478E-01 

0 . 7472E*02 

0“.  1333E+Q0 

0.6726E+02 

0 . 9397E-01 

0.6542E+02 

0.3470E*Q0 

a 

0.4000E+02 

0 . 1 OOOE  +  03 

0 . 4734E+02 
0.2048E-01 

0 . 7732E+02 
0.1810E400 

0.6965E402 

0.1221E400 

0.8693E*02 

0.4429E*00 

9 

0 .4500E402 

0 . 1000E+03 

0.5123E*02 

0 .  2799E-0 1 

0 . 7957E+Q2 
0.2373E+00 

0.7190E*02 

0.1552E+00 

0 .6823E+02 
0.5485E*00 

10 

0 . 5000E+02 

0 . 1000E  +  03 

0 . 5491E+02 

0 . 3774E-0 1 

0.8155E402 

0.3024E400 

0 . 7402E+02 
0.1936E+00 

0.8937E*02 

0.6634E*00 

11 

0 . 5500E+02 

0 . 1000E*03 

0 .5836E+02 
0.5021E-01 

0 . 8329E  +  02 

0 . 3764E  400 

0.7601E402 

0.23S0E*00 

0.9037E*02 

0.7871E*00 

12 

0 .6000E402 

0 . 1000E  +  03 

0.6159E402 

0 .6592E-01 

0.8483E+02 

0 .4593E  +  00 

0.7787E402 

0.2866E400 

0.9126E*02 

0.9191E*00 

13 

0.6500E+02 

0. 1000E  +  03 

0.6460E402 

0.8540E-01 

0.8621E+02 

0.5509E400 

0 . 796 1E+02  . 
0.3459E400 

0.9205E*02 
0 . 1G59E  +  01 

14 

0 . 7000E+02 

0.1000E+03 

0.6740E402 

0 . 1092E+00 

0.8744E+02 

0 .651 1E*00 

0 .8122E+02 

0 .4102E+00 

0.9277E*02 
0  •  1206E*0  Z 

15 

0.7500E402 

0.1000E+03 

0.70008*02 

0. 1379E*00 

0.8855E402 

0.7597E400 

0 . 8272E+02 

0 . 48 1 8E+00 

0.9340E+02 

0 . 1360E+0 1 

16 

0.6000E402 

0 . 1000E+03 

0 , 724IE+02 
0.1720E400 

0.8955E402 

0 . 8765E+00 

0.8410E402 

0.5610E+00 

0.9398E*02 

0.1521E+01 

17 

0.6500E402 

0 . 1000E  +  03 

0.7463E+02 

0 . 2120E+00 

0 . 9046E+02 

0 . 1001E  +  01 

0.8538E*02 

0.6477E*00 

0.9450E*02 

0.1688E+01 
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id 

0.9000E+02 

0. 1000E+03 

0.7668E+02 

0.2583E+00 

19 

0.9500E*02 

0 . 1000E  +  03 

0.7858E+02 
0.31 14E+00 

20 

0. 1000E+03 

0 . 1000E+03 

0.8032E*02 

0.3716E+00 

21 

0.1050E+03 

0 . 1000E+03 

0.8193E+02 

0.4392E+00 

22 

0. 1 100E+03 

0. 1000E+03 

0.8341E+02 

0.5144E+00 

23  * 

0.1150E*03 

0 . 1000E+03 

0.8477E+02 

0.5973E+00 

24 

0.1200E+03 

0. IOOOE+03 

0.8602E+02 

0 .6880E+00 

25 

0. 1250E+03 

0 . 1322E+03 

0.8776E+02 

0.787SE+00 

26 

0. 1300E+03 

0 . 1322E*03 

0.9300E+02 

0.9159E+00 

27 

0 . 1350E+03 

0 . 1322E+03 

0.9916E+02 
0.11  10E+01 

28 

0 . 1400E*03 

0.1322E+03 

0 .  1048E+03 
0. 1403E+01 

29 

0 . 1450E+03 

0 . 1322E*03 

0. 1096E+03 
0. 1820E+01 

30 

0.1500E»03 

0 . 132ZE+03 

0.1138E+03 

0.2385E+01 

31 

0 . 1550E+03 

0 . 1322E+03 

0 . 1170E*03 
0.3114E+01 

32 

0. 1600E+03 

0 . 1322E+03 

0. 1 198E+03 
0.4013E+01 

33 

0 . 1650E+03 

0 . 1322E+03 

0. 1221E»03 
0.50&4E*01 

34 

0.1700E-*03 

0. 1322E*03 

0.1240E*03 

0.6320E+01 

35 

0. 1750E+03 

0. 1322E+03 

0. 1255E+03 
0.7711E+01 

36 

0. 1800E+03 

0 . 1322E+03 

0. 1268E+03 
0.9243E+01 

37 

0. 1650E+03 

0 .  1322E+03 

0. 1278E+03 
0. 1090E+02 

0.9128E+02 

0.1l34E*0l 

0.8657E+02 

0.7421E+00 

0.9497E+02 
0.1W1E  +  01 

0 . 9202E+02 
0.1273E*01 

0.8766E+02 

0.8443E*00 

0 . 9540E+02 
0.2040E+01 

0 . 9270E+02 

0. 1420E  +  01 

0.8866E+02 

0.9541E+00 

0.9579E*02 

0.2223E+01 

0.9332E+02 

0. 1574E  +  01 

0.8959E+02 

0 . 1071E+01 

0.96-15E*02 

0 .241 1E+01 

0 . 9388E+02 

0. 1733E+01 

0 ,9044E*02 

0 . 1 196E>01 

0.9847E+02 

0.2604E+01 

0 . 9439E+02 

0. 1899E+01 

0 . 9122E+02 

0. 1328E+01 

0.9677E+02 

0.2801E+01 

0.9486E  +  02 
0.2071E+01 

0.9194E*02 

0. 1487E+01 

0 .9704E+02 
Q.3002E+01 

0. 1322E*03 
Q.4390E+01 

0 . 1322E+03 

0 . 3787E*0l 

0.1322E+03 

0 .53Z2E+01 

0 . 1322E+03 
0.6710E+01 

0 . 1322E»03 
0.6107E+01 

0. 1322E+03 
0.7641E*0l 

0 . 1322E+03 
0.9030E*01 

0 . 1322E+03 

0 .8426E+01 

0.1322E+0V- 

0.9961E+01 

0 . 1322E+03 
0.11 35E+02 

0. 1322E+03 

0 . 1075E+02 

0.1322E*03 
0. 1228E*02 

0.1322E+03 

0.1367E+02 

0 . 1322E+03 

0 . 1307E+02 

0.1322E+03 
0. 1480E+02 

0 . 1322E+03 

0 . 1599E+02 

0 . 1322E+03 

0 . 1S39E+02 

0.1322E+03 

0 . 1692E+02 

0 . 1322E+03 

0 . 1831E*02 

0 . 1322E+03 

0 . 1771E+02 

0 . 1322E+03 
0. 1924E*02 

0 . 1322E+03 

0 . 2063E+02 

0. 1322E+03 

0 . 2002E*02 

0. 1322E+03 
0.2156E+02 

0 . 1322E+03 
0.2295E+02 

0 . 1322E+03 
0.2234E+02 

0. 1322E+03 

0 . 2388E  +  02 

0 . 1322E+03 

0 . 2527E+02 

0 . 1322E+03 

0 . 2466E+02 

0. 1322E+03 
0 . 2620E+02 

0 . 1322E*03 
0 . 2759E+02 

0 . 1322E  +  03 
0.2698E+02 

0. 1322E+03 
0 . 2852E+02 

0 . 1322E+03 
0 . 2991E+02 

0 . 1322E+03 
0.2930E+02 

0 . 1322E+03 
0.3084E+02 

0.1322E+03 

0.3223E+02 

0.1322E+03 

0.3162E+02 

0. 1322E+03 
0.3316E+02 

38 

0 . 1900E+03 

0.2000E+02 

0.1284E403 

0 . 1266E+02 

0 . 10 IQE+03 

0 . 3294E  +  02 

0 , 9655E *02 

0 . 3207E* 02 

0 . 76  355*02 
0.33395+02 

39 

0 . 1950E+03 

0 . 2000E  +  02 

0 . 1 262E+03 

0 . 1433E+02 

0 . 8690E  +  02 

0 . 3309E+02 

0 , 8685E  +  02 

0 .  3220E *02 

0.61895*02 
0 . 3341 E+02 

40 

0 . 2Q0QE+03 

0 . 2000E402 

0 . 1 2 14E+03 

0 . 1562E+02 

0 . 7742E+02 

0 . 331 7E+Q2 

0 . 8074E*  02 

0 , 3228E*02 

0.5431E+02 
. 0 . 3343E+02 

41 

0.2050E403 

0  *  2000E  +  02 

0 . 1 157E+03 
0.1651E+02 

0 . 7008E+02 

0 . 33205*02 

0 . 7606E  +  02 
0.3234E*02 

0 .4932E  +  02 
0.3344E+02 

42 

0 .  2100E  +  03 

0 . 2000E *02 

0 . 1 097E+03 

0. 1711E+02 

0 . 641  IE  *02 

0 . 3323E*02 

0 . 7204E  +  02 

0 . 32385*02 

0.4559E+02 
0  *  33445+02 

43 

0 . 2 150E+03 

0 . 2000E+02 

0 . 1037E+03 

0 . 1750E+02 

0 . 5914E*02 

0 . 3325E*02 

0 . 6637E*02 

0 . 3241E*02 

0.42625*02 
0 . 33455+02 

44 

0 . 2200E  +  03 

0 . 2Q00E+02 

0.9787E*02 

0 . 1 776E+02 

0.5493E*02 

0 . 3326E*02 

0 . 6492E*02 

0 . 3243E*02 

0.40155*02 
0 . 33455*02 

45 

0 . 2250E+03 

0 . 2000E  +  02 

0 . 9228E+02 

0 . 1793E+02 

0 . 51 32E*02 

0 . 3327E*02 

0 .6 167E*02 
0.3245E*02 

0 . 38055+02 
0.33455+02 

46 

0 . 2300E+03 

0.2000E+02 

0.8697E+02 

0 . 1805E+02 

0 . 481 9E*02 

0 . 3328E*02 

0 » 5859E+02 

0 . 3247E*02 

0.36245+02 

0.33465+02, 

47 

0 . 2350E+03 

0 . 2000E  +  02 

0.8194E+02 

0 . 181 3E*02 

0.4545E*02 

0 . 3328E*02 

0.5569E*02 

0 , 3248E*02 

0.34675*02 
0 . 3346E+02 

48 

0 . 2400E+03 

0 .20005*02 

0.7722E+02 

0.1819E*02 

0 . 4304E*02 

0 . 3329E*02 

0.5296E*02 

0.3249E*02 

0  *  33285*02 
0.3346E+02 

49 

0.2450E+03 

0 , 2000E*02 

0 . 7279E+02 
0.1824E+02 

0.4091E*02 

0.33295*02 

0.5041E*02 

0 . 32505*02 

0.32045*02 

0.33475+02 

50 

0 . 2500E+03 

0 . 2000E+02 

0  *6866E*02 

0. 1827E*02 

0 , 39005*02 

0 . 33305*02 

0.4803E*02 

0 . 325 1E*02 

0.30955*02 
0 . 33475+02 

51 

0 . 2550E+03 

0 .  2000E*02 

0 , 6481E+02 
0.1829E+02 

0 . 37305*02 

0 . 33305*02 

0 . 4582E+02 

0 . 3251E+02 

0 .  2997E  +  02 
0 . 3347E+02 
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APPENDIX  C  ..  LITERATURE  SEARCH  -  BONDING  OF  RUBBER  TO  NONMETALLIC  FIBERS 


The  attached  bibliography  Is  a  result  of  a  literature  search  on  adhesion  of 
nonmetall 1c  fiber  materials  to  rubber.  The  reports  covered  aramld,  nylon, 
carbon,  cellulose,  and  glass  fibers.  In  some  cases,  adhesive  was  applied 
directly  to  the  fibers,  In  other  cases,  bonding  agents  were  added  to  the 
rubber,  and  In  other  cases,  no  adhesives  or  bonding  agents  were  used.  The 
fiber  diameters  were  on  the  order  of  lOum,  which  Is  much  smaller  than  the 
steel  wire  diameter  of  0.032  Inches  (0.8  mm)  which  would  be  used  in  our 
baseline  materials. 

In  reference  C3  (from  the  bibliography),  It  Is  shown  that  bond  strength 
depends  on  the  viscoelastic  properties  of  the  rubber.  Rate  and  temperature 
effects  are  Important.  Adhesion  energy  as  measured  by  peel  tests  for  a  range 
of  rates  and  temperatures  showed  a  reasonable  fit  to  a  "universal"  form  of  the 
WLF  rate-temperature  equivalence  principle: 


lo?Q  [a(T) ]  =  -17 . 4( T-Tq) ,  Tq  =  -90°C 
51.6+T-T 

g 

This  form  Is  compared  to  a  previous  form  (Equation  1)  In  Figures  Cl  and  C2. 

The  difference  Is  very  slight. 

A  direct  comparison  of  the  effect  on  rubber  of  chopped  fibers  of  aramld, 
nylon,  carbon,  cellulose,  and  glass,  with  and  without  bonding  agents,  appears 
In  reference  C6.  The  materials  were  milled  to  get  good  uniaxial  fiber 
alignment,  and  tensile  strength  and  modulus  In  the  fiber  and  transverse 
directions  were  measured.  Samples  with  the  bonding  agent  HRH  showed  the 
largest  Increases  In  modulus  and  strength,  Indicating  a  superior  bond  over 
that  given  by  the  other  bonding  agents  tested.  When  this  bonding  agent  was 
used,  the  strongest  materials  were  those  with  the  aramld  and  nylon  fibers. 
Carbon  and  cellulose  formed  a  lower  class,  with  glass  fiber  reinforcement  the 
least  effective. 

According  to  Reference  Cl,  the  problem  with  the  glass  fibers  may  be  that  they 
are  too  brittle,  and  break  up  In  processing.  Their  studies  showed  that 
chopped  fibers  are  most  effective  when  their  L/D  Is  above  100,  and  that  their 
glass  fibers  after  processing  had  L/D  of  25  to  50. 

The  effect  of  Monsanto  SantowebR  treated  cellulose  fiber  on  cut  growth, 
cutting,  and  chipping  was  examined  In  reference  C7.  Fiber  loadings  of  2  phr 
showed  significant  reduction  In  cut  growth  rate  and  Improvements  In  cut  and 
chip  resistance.  The  reference  does  not  discuss  fiber  alignment.  Processing 
was  done  using  a  Banbury  mixer,  with  no  milling,  so  fiber  orientation  may  have 
been  random. 

A  direct  comparison  of  fiber  pullout  force  In  rubber,  between  Kevlar  aramld 
cords  treated  with  an  epoxy/RFL  adhesive  and  brass  coated  steel  tire  cords  of 
similar  diameter,  was  presented  In  reference  CIO.  It  Is  not  known  If  the 
adhesive  treatment  restricts  the  fiber-rubber  Interface  to  the  outer  surface 
of  the  Kevlar  cord,  or  If  the  larger  total  surface  area  of  all  of  the  Kevlar 
filaments  plays  a  role.  In  any  case,  the  fiber  pullout  forces  for  the  Kevlar 
and  the  steel  cords  were  roughly  comparable,  staying  closer  than  a  factor  of  2 
for  varying  rubber  formulations  and  heat  aging  programs. 
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RELATIVE  TEAR  GROWTH  RATE  R/R(0) 


WLF  FORMULATION 

RELATIVE  TEAR  GROWTH  RATE  IN  SBR  AS  FUNCTION  OF  TEMPERATURE 

LOG10(A(T))  =  -8.86*(T-TS)/ (101.6+T-TS),  TS=TG+20C,  TG=-63C 


Figure  Cl  -  WLF  Formulation 


RELATIVE  TEAR  GROWTH  RATE  R/R(0) 


RELATIVE  DEBOND  GROWTH  RATE  AS  FUNCTION  OF  TEMPERATURE 
"UNIVERSAL"  FORM  OF  WLF  FORMULATION  FOR  RUBBER 

LOG10(A(T))  =  -17.4*(T-TG)/(51.6+T-TG),  TG  =  -90C 


Figure  C2  -  l,Universal',Form  of  WLF  Formulation  for  Rubber 
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